Periodic-orbit theory of universal level correlations in quantum chaos 
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Using Gutzwiller's semiclassical periodic-orbit theory we demonstrate universal behavior of the 
two-point correlator of the density of levels for quantum systems whose classical limit is fully chaotic. 
We go beyond previous work in establishing the full correlator such that its Fourier transform, the 
^\ ' spectral form factor, is determined for all times, below and above the Heisenberg time. We cover 

, dynamics with and without time reversal invariance (from the orthogonal and unitary symmetry 

' classes). A key step in our reasoning is to sum the periodic-orbit expansion in terms of a matrix 

' integral, like the one known from the sigma model of random-matrix theory. 
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Classical chaos is characterized by sensitive dependence on initial conditions, and by the fact that long trajectories 
uniformly fill the available space. These classical features have profound consequences in quantum mechanics: Many 
^ quantum properties of chaotic systems become universal in the semiclassical limit, i.e., they no longer depend on the 
system in question but only on symmetries. For instance, as conjectured in [l| the statistics of energy eigenvalues is 
universal: Neighboring levels tend to repel each other, and the correlation function of the level density p{E), 
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I. INTRODUCTION 



(where the brackets and overbars denote averages over the center energy E) has the same form for all dynamics belong- 
ing to the same symmetry class. The most important classes are those of systems without any symmetries (unitary 
• ■ class; complex Hermitian Hamiltonians) and systems whose sole symmetry is time-reversal invariance (orthogonal 
^5 [ class; real symmetric Hamiltonians). 

If one accepts universal spectral correlations as a fact, a phenomenological prediction for i?(e) can be obtained by 
averaging over ensembles of systems sharing the same symmetries. Modelling the corresponding Hamiltonians through 
matrices, one arrives at random-matrix theory (RMT), as initially proposed by Wigner and Dyson in the fifties 0] in 
the context of atomic nuclei. RMT has since found broad applications in many fields of physics. For the correlation 
function i?(e), random-matrix averages yield the prediction ^, i^, i5|] 



' D/- ^ / ^'''(*^)^ unitary 



^^^^ ^ I -s(e)2 + s'(e) (iSi(e) - isgn(e)) orthogonal 

with s(e) = and Si(e) = /J de's(e'). It is convenient to access R{e) through a complex correlation function C(e) 
which is analytic in the energy upper half plane and connected with the real correlator as i?(e) = ^^^q Re C(e). 

The correlator C (e) has the nice property of being retrievable by Borel summation Q from its asymptotic expansion 



,1—0 \^/ 



(3) 



n=2 

The foregoing power series contains both non-oscillatory terms and oscillatory terms proportional to e^"^; the latter 
are responsible for the singularity (discontinuity at r = 1 of the first or third derivative in the unitary and orthogonal 
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FIG. 1: Bunch of 72 (pseudo-) orbits differing in two 3-encounters and one 2-encounter (a pseudo-orbit is a set of several 
disjoined orbits; see below). Different orbits not resolved except in blowups of encounters 




FIG. 2: Sieber-Richter pair 



case, respectively) in the Fourier transform of the spectral correlator K{t) — ^ (iei?(e)e^'^'^ (the spectral form 
factor). For the unitary class we have C2 = — 5, (i2 = 5 and all other coefficients vanish, whereas the orthogonal class 
involves C2 = — 1, c„ = *'"~^2i^"~^'' ^"-"^ n > 3, d2 = ~ and — (.■n—3y{n-3) n> A. 

This phenomenological RMT result leaves open the question why chaotic systems behave universally, and that 
question we want to address in the present paper. Taking a semiclassical approach one can follow Gutzwiller Q 
and express the level density as a sum over contributions of classical periodic orbits. The correlation function then 
turns into a sum over pairs of orbits. Systematic contributions to that double sum are due to pairs of orbits whose 
actions are sufficiently close for the associated quantum amplitudes to interfere constructively. The task of finding 
correlations in quantum spectra is translated into a classical one, namely to understand the correlations between 
actions of periodic orbits @ . 

The first step in this direction was taken by Berry Q who showed that pairs of identical or mutually time-reversed 
orbits explain the leading coefficient C2 . Higher-order contributions c„ with n > 3 are due to a still not widely known 
feature of chaos, a certain "bunching" of periodic orbits A long orbit has many close "self-encounters" (see 

Fig. [T|) where it comes close to itself in phase space (possibly up to time reversal, for dynamics with time reversal 
invariance). Such an orbit is just one of a closely packed "bunch". All orbits in a bunch are nearly identical, except 
that the orbit pieces inside the encounters are "switched" . 

Based on insight from disordered systems [ll|, the study of bunches of orbits differing in encounters was pioneered 
by Sieber and Richter [T^ who derived the next-to-leading coefficient C3 for time-reversal invariant systems from orbit 
pairs differing in just one encounter where two stretches of an orbit are close (see Fig. [5]). Full agreement with all 
coefficients c„ from RMT was established by the present authors in [l^ . 

In none of these works oscillatory contributions could be obtained (Note however courageous forays by Keating [l^ 
and Bogomolny and Keating [13)) due to the fact that Gutzwiller's formula for the level density is divergent. To 
enforce convergence, one needs to allow for complex energies with imaginary parts large compared to the mean level 
spacing, Ime^ 1. Oscillatory terms proportional to e^"^ then become exponentially small and cannot be resolved 
within the conventional semiclassical approach. 

In [isj we proposed a way around this difficulty that we here want to elaborate in detail. The key idea is to 
represent the correlation function through derivatives of a generating function involving spectral determinants. Two 
such representations are available and entail different semiclassical periodic-orbit expansions. One of them recovers the 
non-oscillatory part of the asymptotic expansion ([3]), essentially in equivalence to [13) [3 [13]; the other representation 
breaks new ground by giving the oscillatory part of ([3]). 

The full random-matrix result also, and in fact most naturally, arises within an alternative semiclassical approxima- 
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tion scheme proposed by Berry and Keating in [19| . That scheme constrains the semiclassical periodic-orbit expansion 
of the spectral determinants det(iJ — H) to be real and to converge for real energy argument. Inserted into the gener- 
ating function the resulting "Riemann-Siegel lookalike formula" for det(£' — H) was shown in to simultaneously 
produce both the non-oscillatory and oscillatory parts of the expansion ^ . Interestingly, the Riemann-Siegel lookalike 
formula parallels an exact composition of the random-matrix averaged generating function originating from a Weyl 
group symmetry [2ll |. 

We want to stress that obtaining both the non-oscillatory and the oscillatory contributions to R{e) also implies a 
derivation of level repulsion: By summing up the asymptotic series ([3]) we access the full R{e) including its behavior 
for small e, i?(e) oc lel*^ (with (3 — 2 ior the unitary and /3 = 1 for the orthogonal class), indicating a suppression of 
small energy differences. To get access to small e the oscillatory contributions (varying on the scale of the mean level 
spacing) are of crucial importance. 

For our semiclassical approach to spectral statistics we have drawn inspiration from a field-theoretical implementa- 
tion of RMT, the nonlinear sigma model. There, the ensemble averaged generating function is ultimately transformed 
into an integral over matrices. In a stationary-phase approximation to that integral, two different saddles are found to 
respectively contribute non-oscillatory and oscillatory terms. The relation between these two contributions bears close 
resemblance to semiclassics. Moreover the topologically different structures of orbit bunches in semiclassics turn out 
to be analogous to Feynman diagrams in a perturbative treatment of the sigma model, with encounters corresponding 
to vertices. That analogy helps to do the sums required in semiclassics and to compact the semiclassical series to a 
matrix integral like the one in the sigma model. This derivation has to be contrasted with a previous field-theoretical 
approach to spectral universality, the ballistic sigma model [1^ . 

The present paper is organized as follows. We start by explaining the generating- function approach (Section|TT]), and 
give the pertinent generalization of Berry's diagonal approximation (Section |TTT|. We then evaluate the contributions 
of bunches of periodic orbits to the generating function and the correlation function (Section lIVp and sum them up 
(Section |V|. Our reasoning culminates in the semiclassical construction of a sigma model lookalike (Section |VI|). For 
pedagogical reasons. Sections IIIIII VII are devoted to the unitary symmetry class, the simplest case. The generalizations 
necessary for time- reversal invariant dynamics are listed in Section [VIII 

The Appendix |^ contains a derivation of contraction rules needed for the semiclassical construction of the matrix 
integral. Further appendices contain ([B| a brief introduction to the supersymmetric sigma model, ([Cl |D|) two replica 
versions of the sigma model and their relation to the semiclassical theory, ([E]) a characterization of orbit topologies 
in terms of permutations, and ([F| the energy average employed in the definition of the correlator and the connection 
between the real and complex correlation functions. 



II. GENERATING FUNCTION 



A. Definition 



All approaches to resolve non-oscillatory and oscillatory contributions start with representing the correlation func- 
tion through derivatives of a generating function. To get there we express the complex correlator C(e) in terms of 
Green functions {E — H)^^ of the Hamiltonian H, 

^ + 1^ - '^)" " - 1^ - - 5 • 

where the superscript -I- denotes a positive imaginary part iry. The average will be understood as over the real center 
energy with an averaging window large compared to the mean level spacing (see Appendix |F|) . The first Green 
function in ([3]) is a retarded one due to the positive imaginary part of the energy argument while the second is an 
advanced one. 

A suitable generating function involves four spectral determinants A(£') — det(i? — H), 

Z(tA, en, ec, eo) ( — ) { ) } ■ (5) 



The four parameters eA,iB,ic,^D denote offsets from the center energy E, made dimensionless by referral to the 
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mean level spacing they are all taken to include a positive imaginary part ir] ^. It is worth noting the symmetry 

Z{eA,eB,ec, cd) = Z{eA, es, -(d, -ec) (6) 

which in the framework of the sigma- model treatment of RMT reflects a Weyl symmetry [2l|. The normalization 
Z(eA, es, CA, es) = 1 also deserves being mentioned. 

The correlation function can be accessed from Z in two different ways. In the standard "columnwise" representation, 
we take two derivatives w.r.t. and es, and subsequently identify all energy increments, 
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Here we used the identity det — exptr In. The notation (||) signals that when equating 6^,63,6(7 and and 
accounting for the signs in ([5]), we have identified the four energy arguments in (O in a "columnwise" way. 

An alternative representation of the correlation function can be obtained if we identify the energy arguments in a 
"crosswise" way (in view of the signs in ([5])), 
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(X) 2(7r/3)- \ • Z7rp ■ ■ zirp A{E 

(x) : = es = e+ , ec = ££) = -e" = -e + i?7 

In ([9|), the ratio of determinants^ becomes unity only in the limit 77 ^ 0, i.e., one can access i?(e) as 
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lim ( 



d^z 



deAdei 



(X) 



(9) 
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B. Periodic-orbit representation 

We are now prepared to semiclassically approximate the generating function in terms of classical periodic orbits. 
Our starting point is Gutzwiller's formula for the trace of the Green function, 

tr(i?+ - H)-^ ^ ~tnp(E+) 1 ^ T,F,e'^»(^^)/^ . (11) 

a 

Herein, the first term involves the average level density p given by Weyl's law a,s p — ^^^.^^ (for systems with two 
degrees of freedom), with the volume of the energy shell. The orbit a is represented by its period T^, action Sa, 
and stability amplitude Fa (defined to include the Maslov phase factor but not the period) 



^ Due to the imposed positivity of the imaginary parts of the energy offsets the determinants involving and ec are associated with 
retarded Green functions. The minus signs in front of eg and in JSJl let the respective spectral determinants be associated with 
advanced Green functions. Note that the present sign convention differs from the one in our letter [iSll . 

^ The above ratio averages to Z(e+, e"*", — e~, — e~) which can also be written as Z(e+, e+, — e~, — e~) = (exp[27riAAf,,]) with AA'',; = 
Nn{E + ~ ^ri{E — 2^)- Here A'^^ denotes the level staircase smeared out over a range rj, Nrj{E) = ^f,Sn{E — E^) with 

9r]{x) = 7r~^Imln(— x + irf) the similarly smeared step function. The function Z{e'^,e'^, —€~, —e~) could be used as a starting point 
for calculating the oscillatory part of the correlator C(e), and then some similarity with the courageous foray of Bogomolny and 
Keating [TtII mentioned in the Introduction would arise. Indeed, replacing A'^,, by the sum of its energy average and a fluctuating part, 
Nri = (Nr,) + N^'^"'^, and treating N^"-"'^ in a Gaussian approximation we get (exp[27riAAf,,]) ~ e^'^ expl-2n^ {(AN^"-"'^)'^)]. The second 
factor contains the number variance in the exponent; the first factor signals that the oscillatory part of the two-point correlator becomes 
accessible. This rough argument anticipates the diagonal approximation below. 

^ We disregard orbits which are multiple repetitions of a shorter orbit since these are exponentially suppressed in the limit of long periods. 
As regards orbits incorporating multiple approximate traversals of short orbits, the respective corrections were shown to be immaterial 
for the non-oscillatory part of the spectral correlator (see and can similarly be shown to cancel in our four-determinant generating 
function Z. Interestingly, such orbits do contribute to the averaged product of two spectral determinants [Tsll . 
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To proceed to a semiclassical approximation of the spectral determinant, we integrate and exponentiate in pT|) as 
A(£'+) cx exp ^ dE'tT-^TZTfj^ ■ The resuhing exponent inherits the smoothed number N{E) of eigenvalues below 
E and an orbit sum from the trace formula (jlip . 

A{E+) oc exp ( - i7rNiE+) - ^ ^^.e'^"'^^)/'') = e-''^^^^^) ^ i^^(_i)"AeiS^(£+)A . (12) 

a A 

The last member of the foregoing chain of equations has the exponentiated orbit sum expanded into a sum over 
pseudo- orbits, i. e., non-ordered sets of periodic orbits A [l^. In that pseudo-orbit sum, tia is the number of orbits 
inside A, Sa the sum of their actions, and Fa the product of their stability amplitudes Fa. The sum over A includes 
the empty set, whose contribution is unity. (Note that in the foregoing expressions we suppressed a proportionality 
constant due to the lower limit of the iJ'-integral. That factor is irrelevant since it cancels in the generating function. 

Eq. (jl2p applies to energies with a positive imaginary part. For energies with a negative imaginary part complex 
conjugation yields 

A{E-) = A{E+y (X exp [iT:N{E') - ^ _p*g-iS„(£;-)/a^ 

a 

A 

For inverse spectral determinants, the minus sign in the periodic-orbit sum and hence the factor (—1)"-* in the 
pseudo-orbit sum disappear, 



A{E+)-' = (A(£;-)-i)* cxexp(i7riV(£;+)+^Fae' 



s„(£;+)/ri 



A 

The four spectral determinants in Z now yield semiclassical asymptotics as a fourfold sum over pseudo-orbits. 



^g-i7rlV[E-eB/(27rp)] ^ p^^-iS b[E -es / i27rp)]/h 
B 

^g--i7rlV[£;+£c/(27rp)] ^ ^_^-^nc ^iSclE+ec / {2TTp)]/h 
C 

^^i7T'NlE-en/{2Tvp)]^p*^f^_^yiDQ-iSo[E-eD/{2TTp)]/h\ ^ ^-j^g-j 
D I 

where all offset variables tA^^Bi^CT^D have positive imaginary parts 77. Here the superscript 1 signals that (jlSp is a 
semiclassical approximation based on the Gutzwiller trace formula which requires 

7?>1 (16) 

for convergence such that terms vanishing in that limit (such as terms of the order e"'') will be lost. The complement 
missed by Z^^^ becomes important for real energies (a manifestation of the Stokes phenomenon of asymptotic analysis) 
and will be introduced below. 

Upon expanding as N{E ± ^l^) ~ N{E) ± ^ and S{E ±^)/n^ S{E) ± ^ = S{E) ± Te/Tu where Th = 

2'kTip — ipj- is the Heisenberg time, we simplify the quantity Z*^^^ to 



^(1) ^ ^i(e^+eB-ec-ei,)/2 X ^ ^ F^F* Fci^^ (-1)"^+"" 



A,B,C,D 

y^Ql(SA[E)-SB{E) + Sc{E)-SD{E))/h^l(TAtA+TBtB+Tcf-C+TDtD)/TH\j _ (27) 

The foregoing sum over pseudo-orbit quadruplets combines classical chaos with quantum interference due to the 
phase factor e"^'^/'' = Q^iSA{E)-SBiE)+SciE)-SDiE))/h ^ -p^j. jj^q^^ quadruplets the interference is destructive since for 
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?i — > the phase factor e'^"^/^ oscihates rapidly as the energy varies and therefore vanishes after averaging. For 
constructive interference, the action mismatch AS" must be small, at most of the order of h; to within such a quantum 
mismatch the cumulative action of the pseudo-orbits BUD must coincide with the cumulative action oi AuC. Such 
action matching happens systematically only if the orbits in B U D either coincide with those in AU C (and thus 
constitute the "diagonal" contributions, see Sect. IIIip or differ from them only by their connections inside encounters 
("off-diagonal" contributions) as in Figs. Il|2l Indeed, then, the phenomenon of orbit bunching is crucial for spectral 
universality. 

Now imagine Z^^^ calculated from P7)) and then analytically continued to real energy offsets, ry — > 0+. If instead 
of we had the full generating function we could obtain both the non-oscillatory and the oscillatory contributions 
to C(e) by taking derivatives w.r.t. and cb and identifying the energy increments either in the columnwise way 
([7]) or the crosswise one ©■ However, for Z^^^^ the columnwise procedure ([7]) cannot yield oscillatory terms since it 
makes the (Weyl) exponential e''^'^'^+'^^~'^'^~'^")/^ collapse to unity. In contrast, with the crosswise representation Q 
the signs of ec and en are flipped and the Weyl factor turns into e^"^. Hence the crosswise representation recovers 
only the oscillatory Fourier component proportional to e^'^ . 

The random-matrix result suggests that to get the full correlation function one has to add both results 



(II) deAdes 



i?(e) = -Re 2— ^ + 2— ^ _ + - . (18) 



(X) 2 



Here the first and the third summand are due to ([7]) while the second summand is due to (fTO|) ^. Of course it is vital 
to justify this additivity independently from RMT, a task to be attacked below. 

7(1) I ^2 7(1) I 

The "crosswise" term q^^q^^ | ^^-j differs from the "columnwise" term -§^^q^\(j^ because ec and e/j are interchanged 
and flipped in sign. This interchange and sign flip may as well be performed on the level of the generating function, 
i.e., we can write 



deAdes 



(19) 



(II) 



Z^^\eA,eB,ec,eD) = Z^^\eA,eB,-eD,-ec) (20) 

Instead of adding the crosswise representation in we can use Z^^^ -I- Z'^2) ^s a generating function, and get 
i?(e) in full through the columnwise (or crosswise) representation alone. Therefore no extra work is needed for the 
oscillatory part of the two-point correlation function, once the "standard" part Z^^^ of the asymptotics of Z is found. 
The resulting asymptotics of the generating function 

Z-Z(i)+Z(2) (21) 

respects the Weyl symmetry and remains valid when the arguments are real (77 = 0). We note that if one wants to 
evaluate Z'^2) using the Gutzwiller trace formula, the sign flip in ec and eo means that the corresponding arguments 
of Z*-2) should be taken with large negative imaginary parts. 



C. Riemann-Siegel lookalike 

As pointed out by Keating and Miiller in f20'| , the additive structure (PT|) of the generating function most naturally 
arises with the help of the Riemann-Siegel lookalike formula [l^. 

The essential idea is to eliminate a drawback of the semiclassical approximation (|12ll3p of the spectral determinant, 
namely its failure to manifestly preserve the reality of the energy eigenvalues, i. e., the unitarity of the quantum 
dynamics. In the limit r] ^ 0, A(£'+) and A(£'^) must become real and identical to each other. That property 
is not manifest in but may be imposed on Eqs. (|12I13[) . By imposing this property, it was shown in [l9| that a 
duality arises between "long" pseudo-orbits (sum of periods of contributing orbits larger than half the Heisenberg 
time Th) and "short" ones (cumulative period smaller than Th/^); the overall contribution of long orbits is the 



* Critical readers may miss the i from JTOj. This summand does not show up here since it is non-oscillatory and cannot be retrieved 
from the crosswise representation. 
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complex conjugate of the contribution arising from short orbits. The ensuing Riemann-Siegel lookahke expresses the 
real spectral determinant A{E) as a sum over pseudo-orbits with durations shorter than Tii/2, 

A(£;) oc e~"^(^) FA(-l)"^e'^^/'' +C.C.. (22) 

A (Ta<Th/2) 

To evaluate the generating function Z one now uses the Riemann-Siegel formula for the two spectral determi- 
nants in the numerator. For the determinants in the denominator we are not allowed to use Riemann-Siegel as the 
non-zero imaginary parts with appropriate signs in eA,^B are crucial for the definition of Z, and we just stick to 
Eq. (jl4p . Since each determinant in the numerator becomes the sum of two mutually complex conjugate terms, 
Z has altogether four summands which involve the four phases ±i7r |lV(i? -1- ec/CSTrp)) — (i? — e_D/(27rp))] and 
±i7r pV (i? -I- ec/(27rp)) -\- N [E — ejj / {2tt'P))\. The two latter contributions oscillate rapidly as functions of E and 
thus vanish after averaging. 

The two surviving additive terms yield the high-energy asymptotics of the generating function Z ^ Z^^^ + Z'^' 
found in the previous subsection, according to ([T71 [^D|) , except that the sums over C and D are now restricted to 
pseudo-orbits with durations shorter than Tu/'i.. To make that difference irrelevant it is convenient to evaluate the 
pseudo-orbit sum for large imaginary parts 77 of the dimensionless energy offsets in the arguments of the action, 
Sa{E) Sa{E + ir]/2Trp), see ([T5)) . That artifice makes sure that in the relevant integrals over orbit periods to 
be encountered below the integrand falls to exponentially small magnitude as the integration variable rises to near 
half the Heisenberg time. If all four energy offsets €a,b,c,d are thus provided with large positive imaginary parts, 
the two additive pieces obtained through the Riemann-Siegel lookalike become respectively Z^^\eA,£B,£Cj^D) and 
Z(^^(eyi, es, Cp, e^i) = Z^^^€a, es, —^*d, —^c)- Complex conjugation of the last two arguments of Z^"^^ is needed to 
enable its calculation in terms of periodic orbits according to (|17p . At the final stage we use analyticity to let 77 go to 
zero; the sum of Z^^^ and Z'-^'> yields then the above asymptotic representation (PT|) . 

To obtain a quantitative understanding of spectral statistics, we now have to evaluate the quadruple sum over 
pseudo-orbits in Z'-'^' , taking into account orbits in B LI D that either (diagonal contributions) coincide with or 
(off-diagonal contributions) at most differ in encounters from the orbits in AU C. 



III. DIAGONAL APPROXIMATION 



In the diagonal approximation we effectively neglect all correlations between orbits that are not identical. If we 
represent the generating functions as a sum over A, B, C, D as in P?|) this means we should keep only contributions 
where the orbits in ^ U C are repeated inside B U D. However we can save labor if we write the four spectral 
determinants entering our generating function like the second member in (jl2p . i. e., with the help of exponentiated 
orbit sums. We thus obtain the equivalent representation of the generating function 



^(1) ^ei(eA+.B-.c-.z,)/2/ Jl^^ 

\ a / 

wherein each periodic orbit participates with the factor 



(23) 



exp 



— fsi 



(24) 



If we assume that contributions of different periodic orbits are uncorrelated, Z^^'^ factorizes into a product of energy- 
averaged orbit factors Za, 



diag 



n 



(25) 



Inasmuch as the stability coefficient Fa ^ Q-^Ta,/2 -g gjj^a,ll for long orbits we may expand the exponential in Zq 
and limit ourselves to the first nontrivial term, [za) ~ 1 -I- \Fa\^ fAcfsD ~ e'^"' fAcfso- here terms proportional to 
Q±iSiE)/h g^^^ g±2iS(_E)/ft hg^vg been eliminated by the energy average. We so obtain 



diag 



(26) 
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This result allows for an intuitive interpretation: If we assume that the non-identical orbits are uncorrelated the 
exponent in the product of all orbit factors Za, namely <I> = [FaG'^'^''^'^/^ f^fj + F*e~'-^'^^^'>/^ fgjj'\, effectively 
becomes a sum of many independent random contributions; due to the central limit theorem $ must then obey 
Gaussian statistics such that (e*) = e^* which just yields (j26|) . 

To bring our result to a more explicit final form, we need to invoke ergodicity in the form of the sum rule of Hannay 
and Ozorio de Almeida [1, H^l : To sum a smooth function of the period T over an ensemble of orbits weighed with 
IFqP, we may as well compute an integral over periods weighed with 

Y.\Fa\\---)^ rf^(...); (27) 
here Tq is some minimal period starting from which orbits behave ergodically. Employing the sum rule (|27p we get 



Since for T ^ Th the integrand is cx {T /ThY can set the lower integration limit to 0, accepting an error of the 
order {Tq/Th)"^ ■ Then using ^ (e^""^ — e'''^) = In ^ we arrive at the diagonal part of the generating function 



Z\ > = c^(<^A+eB-<ic-en) T ^u) yep, -p ^u) ^ ^^g) 



The diagonal approximation yields an estimate for the correlation function C(e). Taking into account both -^^iag 
and its counterpart -^diag defined in Eq. ([20|) we obtain 

deAdeB '(II) 2"2(ie)2 2(ie)2 ' ^^^^ 

This result already reproduces the full RMT prediction for the unitary symmetry class under discussion. It remains 
to explain, within our semiclassical approach, why all off-diagonal contributions do in fact vanish. 



IV. OFF-DIAGONAL CONTRIBUTIONS 



Off-diagonal contributions arise from quadruplets of pseudo-orbits where not all orbits of ^ U C are simply repeated 
inside BUD; sonic orbits of A U C and BUD differ in encounters. By studying these encounters in detail we shall 
show that for systems without time-reversal invariance the off-diagonal contributions cancel mutually in a nontrivial 
way. This also sets the stage for a generalization to time-reversal invariant dynamics where this cancellation does not 
occur, and off-diagonal quadruplets of pseudo-orbits give rise to terms of higher order in 1/e. 



A. Encounters and orbit bunches 



Long periodic orbits often come close to themselves, and one another, in phase space. We speak of an /-encounter 
whenever I stretches of one or more periodic orbits in a chaotic system approach in this way. Each stretch begins 
and ends at phase-space points to be called "entrance ports" and "exit ports" . These ports can be uniquely defined 
by setting an upper bound for the separation between the encountering stretches. The parts of an orbit in between 
encounter stretches will be called "links" . 

Now the crucial point is that the orbits and pseudo-orbits in hyperbolic motion exist in bunches. The members 
of a bunch differ noticeably only in encounters where the ports are connected differently by encounter stretches. In 
contrast the links are very nearly identical. All members of such a bunch have almost identical actions. If two of 
them are taken a.s AuC and BU D in Eq. ([TT]) their contributions interfere constructively and need to be taken into 
account for a full understanding of spectral statistics. 

The phenomenon of bunching is a consequence of hyperbolicity. For a system with two degrees of freedom^, 
hyperbolicity means that deflections between phase-space points in the three-dimensional energy shell can be split 



In the present paper we shall stick to two-dimensional systems for convenience but the generalization to systems with arbitrarily many 
degrees of freedom along the lines of [13 . 125II is straightforward. 
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into three components: one pointing along the flow, one unstable component u, and one stable component s. As two 
close-by phase-space points move, the unstable component of their separation grows exponentially whereas the stable 
component shrinks exponentially. In the limit of large times the rate of growth or shrinking is given by the Lyapunov 
exponent A as 

u{t) - u(0)e^* , s{t) - s(0)e-^* . (30) 

As a consequence, given any two encounter stretches that are long compared with the characteristic size of the 
system and close enough to linearize the motion as in (|30p we can draw a piece of a physical trajectory that starts 
close to the first stretch and ends close to the second stretch. Its separation from the first stretch must point along 
the unstable direction at the respective phase-space points (such that it approaches the first stretch for large negative 
times) whereas the separation from the second stretch must point along the stable direction (such that it approaches 
that stretch for large positive times). The new piece of trajectory is thus located where the unstable manifold of one 
encounter stretch of A U C intersects the stable manifold of the other stretch [2^ . As a consequence we can reconnect 
the ports of each encounter in any way we like while leaving the links in between the encounter stretches practically 
unchanged. (Note that subsequent to the construction just explained exponentially small corrections to the encounter 
stretches and links are needed to connect them to continuous orbits.) 



B. Structures of pseudo-orbit quadruplets 

The pseudo-orbit quadruplets can have many different topologies some of which are shown in Fig. [3] For their 
classification we introduce the notion of a structure by which we mean a diagram like the ones in Fig. [3] but with the 
encounters and their ports numbered and each periodic orbit assigned to a particular pseudo-orbit. All structures can 
be generated if we proceed as follows: 

SI: Encounters - We consider any number V of encounters. For each of the encounters a ~ 1, 2, . . . , V^, we take into 
account all possible numbers of stretches l{a) = 2,3,4, . . .. The l{a) entrance and exit ports of each encounter 
are numbered from 1 to l{a) in such a way that the i-th entrance port is connected to the i-th exit port before 
reconnection, and to the (i — l)-st exit port after reconnection; the first entrance reconnects to the last exit. 
The encounters can thus be depicted as in Fig. [3] where full and dashed lines denote the encounter stretches 
before and after the reconnection. 

S2: Links - We connect the ports by links in all possible ways. Each link must lead from an exit port to an entrance 
port of either the same or a different encounter. 

S3: Orbits - We include all possible distributions of the orbits among the pseudo-orbits A, B, C, D. The graph 
formed by all links and encounter stretches before reconnection falls into a certain number of disjoint parts each 
standing for a periodic orbit. Each of these can be included in either A or C. After changing connections inside 
the encounters the graph consists of a (possibly different) number of disjoint orbits which have to be distributed 
between B and D. 

Thus equipped we shall split the sum over quadruplets in (fT7|) into a sum over structures and a sum over quadruplets 
associated to each structure. The sum over quadruplets associated to each structure will be done using ergodicity. 
We shall determine the likelihood for having encounters inside and between periodic orbits; this likelihood depends 
on the length of the links and the phase-space separation between the encounter stretches. The sum can then be 
replaced by suitable integrals. The subsequent summation over structures becomes a purely combinatorial problem. 

When doing the sums and integrals mentioned, we have to be aware of a certain overcounting. Since our definition of 
a structure involves an ordering of encounters and ports, each physical quadruplet can be associated with a structure 
in several different ways. To begin with, the encounters of a given quadruplet can be numbered in any one of V\ 
different orders. Then, in each encounter any of its l{a) entrance ports can be chosen as the first one. Note that as 
soon as the first entrance in an encounter is chosen, the numbers of all other encounter ports of the quadruplet become 
uniquely defined by our numbering scheme. In all, we get VI 1{(t) alternatives either leading to different structures 
or to the same structure but with different numberings and thus to a different set of the integration parameters. In 
both cases after summation over structures and integration over the continuous parameters the same quadruplet is 
taken into account several times. In order to compensate this overcounting, the contribution of each structure needs 
to be divided by V\ K'^)- Details can be found in Appendix lEl where a formal definition of structures in terms of 
permutations is given. 
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FIG. 3: Physically distinct diagrams contributing to the leading (L — V — 1) and next-to-leading order {L — V = 2) of the 
■i- expansion of Z^^\ Only the ones with shaded backgrounds arise for the unitary class (only parallel encounters and single 
sense of traversal permitted, as indicated by arrows); they give mutually cancelling contributions in each order, both for the 
unitary and the orthogonal classes. The numbers n and n' of orbits in respectively AuC (full lines) and BU D (dashed lines) 
are written above each column as n : n'. 
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FIG. 4: Topological picture of an encounter where connections 1 — 1', 2 — 2', 3 — 3', 4 — 4' are replaced by 1 — 4', 2 — 1', 3 — 2', 4 — 3'. 




FIG. 5: Simplest (pseudo-)orbits. 



C. Classical properties of orbit bunches 

We now follow [l^, [3| to investigate the classical properties associated to each structure, and then evaluate their 
contribution to the generating function. We start with the example of orbits differing in a 2-encounter (see Fig. [5]). 
First we need coordinates measuring the phase-space separations between the two encounter stretches. We therefore 
place a transversal Poincare surface of section at an arbitrary position inside the encounter and consider the two points 
where the encounter stretches pierce through that section. The separation of the piercings can be decomposed into a 
stable component s and an unstable component u [1^ [13, HI] ■ These coordinates determine the difference between 
the action of the original orbit and the cumulative action of the two partner orbits as AS* ~ su (see also [12, 29]). 
Moreover, if we define the encounter as the region with phase-space separations ]s], \u\ < c (where c is arbitrary but 
classically small) the time until the end of the encounter can be estimated as In and the time since the beginning 

of the encounter is obtained as ^ In -ily . Summation thus yields the overall duration^ of the encounter ionc ^ j In jf^ 
p^ . For the semiclassically relevant encounters we need AS" = su ^ h, i.e., the encounter duration is of the order of 
the Ehrenfest time Te = In Encounters are thus considerably shorter than the orbits, whose periods must be of 
the order of the Heisenberg time Th oc h~^. 

The likelihood of finding encounters in a periodic orbit can be determined using ergodicity, i.e., the fact that long 
trajectories and (ensembles of) long periodic orbits uniformly fill the energy shell [30]. Ergodicity implies that the 
probability for finding two piercings through a Poincare section with separations t, s, u in differential intervals dt, ds, du 
is uniform and reads , with the volume of the energy shell. Statistical independence of the two piercings is 

assumed here; that assumption is an allowable one since for orbits with a period of the order Th the links between 
the encounter stretches typically have durations much longer than the classical timescale needed for the decay of 
correlations. 

Now integration of the probability density l/f2 yields the number density wt{s,u) for finding, inside an orbit of 
period T, 2-encounters with stable and unstable separations s and u. That density was determined in [13 , [3ll ] as 

. X T fdt 

Utcnc{s,U) 

Here the factor T is the period of the orbit (originating from integration over the time of the first piercing). J dt 
is an integral over the time of the second piercing reckonned from the first one, with the following restriction: The 



We note that for general hyperbolic systems (where Eq. I I30I I is only valid in the limit of large times) our formula for tcnc is just an 
approximation. The treatment can be made more rigorous if one avoids this approximation, by following the lines of [25j and Appendix 
B of [Hi. 
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time separation between the two piercings must be at least tone so that in between the piercings the orbit can go 
through the rest of the first encounter stretch, a hnk of positive duration, and the beginning of the second stretch. 
This reasoning leads to J dt ^ T — 2tcnc, but we prefer to just write / dt. 

Now let us consider more complicated structures with, say, V different encounters inside or between, say, n orbits 
of periods Ti,T2, . . . ,Tn. The encounters are labelled by an index a running from 1 to 1^ and the trth encounter 
may have l{a) stretches. Drawing a Poincarc section in each encounter we face altogether L — '^^l{<j) piercings, 
l{a) ones for the crth encounter. The separations between the piercings belonging to one encounter are characterized 
by I {a) — 1 stable and unstable coordinates So-m, Uo-m with to = 1, 2, ... , I (a) — 1 (see [l3 | for the precise definition). 
These coordinates determine the action difference as AS* ^ SamUam and the duration of the cr-th encounter as 

tenc.cr = X — ^ — {maK , \u — Tf H ■ By foUowiug the Same steps as for (pij) the number density of encounters is 
generalized to 

^ftftfc ,32, 

with s = {so-m}i u — {ua-m}- Here the orbit periods arise from integrating over the time of the first piercing in 
every orbit, and the integral / d^~^t goes over the possible times of the remaining piercings reckoned from the first 
one. 



D. Contribution of orbit bunches to Z 



We can now investigate the impact orbit bunches have on spectral statistics. We thus evaluate the generating 
function Z'^) as given by the fourfold pseudo-orbit sum in (fT7|) . taking into account quadruplets with AU C and 
B U D differing in any number of encounters. As before AU C and B U D may also coincide in any number of 
component orbits, and therefore ^^[^g ^^^^^ arises as a factor. We can thus rewrite (fTT]) as 

Z^'^-Z^^la + Z^j^), (33) 

with picking up non-empty pseudo-orbits AU C and BUD that differ in encounters while no longer comprising 
identically repeated orbits, 

^ {FaFbFcF^ {^_lYc+nD^i/\S/n^i(TA<iA+TB<^B+Tcec+TDtD)/TH\j 



A,B,C,D 
ditf. in one. 



E (\Fa\'^\Fc\'^ i^_'^Yc+nDQi^S/h^i{TA€A+TBiB+Tcf-c+TD(-D)/TH\ _ (34-) 



diff. in cnc. 



The off-diagonal part of the generating function can now be written as a sum over structures, and then over 
all quadruplets of pseudo-orbits pertaining to each structure. We thus obtain 

(-1)"^+"° 



4'J= E E 



VIU 1(a) 

structures (A, s,c,D) ixcr v / 

e structure 



x/|i;^^|2|^^|2giAS/?igi(TAeA+TBeB+Tc£c+TB£B)/TH\ (-35^ 



where as explained above we divided out the number of possible orderings of encounters and stretches, V\ Y[„ since 
otherwise each quadruplet would be counted multiple times. The sum over pseudo-orbits permitted by a structure 
can be done by summing over all choices for the n component orbits pi,p2, ■ ■ ■ ,Pn of AU C and integrating over the 
L — V — J2aiK'^) ~ 1) pairs of stable and unstable coordinates characterizing the V encounters, with the number 
density ([5^ as weight, 

^il' = E wnT^ E ni^-i^(/^^-"^^^-"" (36) 

structures i- i-cr \ / ,. ■ ■ ,Pti i—1 
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Each of the n orbit sums can now be done using the sum rule of Hannay and Ozorio de Almeida, Y^p. \Fp-\^ 
Importing the density w from p2p we see the orbit periods in w cancel against those brought in by the sum rule and 
face an L-fold time integral, over the n periods and L — n piercing times. 

To simplify the foregoing expression for we change the L-fold time integral just mentioned to an integral over 
the L link durations. The Jacobian of that transformation is unity. We thus obtain 

(1) _ (_l)nc+»l, (jrr Ma)^l Ma)~l ^"-i /M f37) 

striirtiirrs 1 Icr \ / \ \„— I-' CIlL,o / 



where we inserted compensating powers of the Heisenberg time. The summand coming from a structure in p7p then 
becomes a product over all links times a product over all encounters. It is easy to evaluate the pertinent factors. 
Each link appears precisely once in ^ U C as well as in B U Z?. The integral over its duration picks up the integrand 
1 ei(e^ or c+£B or i3)tii„k/T« g^^j coutributcs 

link factor ; (38) 

CA or C + es or D 

clearly, the index alternatives have to be decided according to which pseudo-orbit the link pertains to before recon- 
nection {A or C) and after reconnection {B or D). 

We now turn to the factor coming from an Z-encounter. In the original pseudo-orbit AU C, each of the I stretches 
may belong to either A or C. We denote by Ia the number of stretches belonging to A and by Ic the number of 
stretches belonging to C, with Ia + Ic = l- When switching connections inside the encounter these stretches are 
replaced by I new stretches differently connecting the same ports as the old stretches. These new stretches belong to 
either B or D, and the corresponding numbers are denoted by Ib and Id with Ib + Id = I- The encounter now yields 
a contribution lAtenc to the duration of pseudo-orbit A; analogously for B, C and D. Hence the encounter draws the 
part e'(''*<='*+'^^s+^<^'^<^+'"'^")*=-^/^« from the phase factor e'(^'i^'*+^seB-KTcec-HTD£D)/rH ^ Collecting everything 
from p7p pertaining to the encounter under consideration we get the dimensionless term 



For the calculation of the remaining integral we refer to [IJ, |32| and here but briefly sketch the argument. The key 
is to Taylor expand the second exponential and show that only the term linear in icnc survives in the limit h —^ 0. 
In that surviving term the encounter duration tone cancels and the factor i^Ia^a + Ib^b + Ic^c + Id£d)/Th can be 
pulled out from the integral. Using Th = (^)'~^ / ^'"^■s d'"^we'^™ s^u^/h ^ (^Jd^ e'"")'"^ ^ 1 we get 

^IIa^a + Ib^b + Ic^c + Id^d) encounter factor . (40) 

The sum over structures (j37p thus simplifies to 



7(1)^ 1 / .yic+nD Y\cnc'^i^A^A + lB<^B + lc^C + Id^d) 



A further simplification can be achieved after inspecting more closely the encounter contributions. For an /-encounter 
of some structure we include the overcounting factor j and face the contribution ji(/^e^ -I- Ib^b + Ic^c + Id^d)- 
Each of the I structures obtained from the given one by cyclic renumbering of the stretches of the selected encounter 
contributes identically, and altogether these / structures contribute i{lA^A + Ib^b + Ic^c + Id^d)- We now claim that 
the same result is obtained if we choose to assign to each structure the encounter factor i(e^ or c + 01 d) where 
the index alternatives must be decided according to the placement of the first entrance port] e.g., for a structure 
where the entrance port of the first encounter stretch belongs to the pseudo-orbit A before reconnection and to 
D afterwards we would write i(eA + eu). Taking a sum over our group of I equivalent structures we then obtain 
X] or c = Ia^a + Ic^c since the first entrance belongs to A and C in I a and Ic members of the group, respectively. 
Similarly we have 'Y^tB or d ~ Ib^b + Id£d- Consequently, the encounter contribution ji(Z^e^ + Ib^b + Ic^C + Id^d) 
can be replaced by i{eA or c + or d) since after summation over the group of equivalent structures associated with 
the encounter the two expressions give the same result. Eq. (j41[) can thus be replaced by the simpler formula in which 
the numerator refers to the first entrance port of each encounter, 

^ ^( lync+np IIcncKeAorC +eBorD) |,^2) 



stru'^urcs ^' nii„ks(-i(eA or C + ^B or d)) 



1 

2 



FIG. 6: (a) 2-encounter with ports labelled, exit ports primed; (b, c) two ways of connecting ports with links. Pre-reconnection 
encounter stretches full lines and post-reconnection ones dashed 



A symmetry between link and encounter factors thus results and entails an interesting property of the foregoing 
"master formula" for Z^^: Apart from a factor —1, each encounter factor is cancelled by the link factor for the link 
ending at the first entrance port of that encounter. Indeed, the first entrance port and the link ending there belong 
to the same pseudo-orbit before the reconnection, either A or C, and likewise after the reconnection, either B or D. 
After that cancellation, the numerator of the e-dependent fraction becomes (—1)^ and only the denominator remains 
e-dependent. That property will be important for proving the cancellation of all encounter contributions in Sect. IVl 
According to (j42|l the contribution of a structure to the generating function is of the order -^ikrv- The sum over 
structures therefore generates a power series in ^. Low-order diagrams can easily be drawn; for instance, Fig.[3]depicts 
all diagrams contributing in the leading two orders, - and both for the unitary and orthogonal class (apart from a 
further diagram containing two copies of Fig. [5] and a diagram containing one copy of Fig. [5] and one Sieber-Richtcr 
pair). 



E. Example: Structures for single 2-encounter 

It is well to illustrate the formula ([1^ for the simple example of a single 2-encounter with two entrance and exit 
ports, depicted in Fig. [6] (a); the full and dashed lines respectively show the encounter stretches before and after 
reconnection. We can connect the ports by links (stage S2) in L! = 2 ways; the two resulting diagrams are shown in 
parts (b) and (c) of Fig. [6] . 

We first deal with Fig. [S] (b) where we face a single pre-reconnection orbit 70 incorporating all ports and links. 
After the reconnection the graph disconnects into two orbits 71 and 72 each of which contains one entrance and one 
exit port as well as one link; the labels 1 and 2 signal the number of the respective entrance port. 

At the stage S3 we generate structures by dividing the orbits among the pseudo-orbits. Since we deal with three 
orbits there are 2^ distinct such allotments; four of these are 

1. A = {7o}, S = {71,72}, C = i? = 0, 

2. yl-{7o}, i5 = {71,72}, C = B = 0, 

3. = {70}, B = {71}, ^ = {72}, C = 0, 

4. A = {70}, D = {72}, B = {71}, C = 0, (43) 

and four similar structures have A interchanged with C. Consider the first structure in the foregoing list. The entrance 
port 1 belongs to 70 € A before, and to 71 G -B after reconnection. Therefore the encounter factor in the numerator 
reads i (e^ -I- cb)- Both links belong to A before and B after reconnection; hence the denominator is [— i [ea + es)]^- 
The factors in are calculated using V — 1, tiq = jt-d = 0; the structure thus yields Zi = ^^^^^ ■ The contribution 
of the second structure differs by the replacement B D; there are two orbits in D and none in C such that the sign 
is unchanged and we get Z2 = . In the third structure the first entrance port belongs to A before and B after 

reconnection. The two links are in A before reconnection; after reconnection one of them goes to _B, another one to 

D. The sign factor is negative since nr — 0, tid = 1. Hence we get Z-x — —-, — 't'^^+^-s) — ^ = ! — Similarly, the 

fourth structure gives Za = —-, — '('^-^+^°) — ^ = ! — . Obviously, Z-i cancels with Zn and Za with Z^ . 

The reader is invited to repeat the calculations for the remaining four structures and again find full cancellation. 
Finally, a patient reader will check that Fig. [5] (c) leads to 8 more structures whose contributions coincide with those 
of Fig. [S] (b) and also cancel pairwise. 
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V. CANCELLATION OF OFF-DIAGONAL CONTRIBUTIONS 



We now show that for systems without time-reversal invariance all off-diagonal contributions mutually cancel. Our 
proof is based on two cancellation mechanisms: Each structure is cancelled either by a structure in which the orbits 
are partitioned differently into pseudo- orbits or in which one link has been removed. 



FIG. 7: Case 1 of cancellation: Structure involving orbit with single encounter stretch 



Different partition: To start we consider structures where one of orbits (say, in A) only contains one encounter 
stretch with two ports and a link, as in Fig. [71 and the stretch is the first in its encounter. The ports and the link then 
also have to belong to the same partner orbit (included, say, in B). Then according to Eq. ([1^ the contribution of 
the structure contains a factor i{eA + cb) from the entrance port of the stretch and a factor from the link; 

these link and stretch factors multiply to —1. 

Now we compare with a structure where the orbit from Fig.[7]is included in C rather than A. Then the above factors 
are replaced by i(ec -I- es) and ^^id the product is again —1. However, there is an extra factor —1 because 

the number of orbits in C is increased by one. Hence the two structures considered have cancelling contributions. 

The same argument goes through if B is replaced by D, or if the orbit in Fig. [7] is a partner orbit; in the latter case 
its inclusion in B and D yields mutually cancelling contributions. In all cases we are free to drop from Eq. (j42p all 
structures where an orbit includes just one link and one encounter stretch which is the first in its encounter. 



!0 



(a) (b) 

FIG. 8: Case 2 of cancellation: Structure involving 2-encounter and Z-encounter cancels against structure with the two encounters 
merged to an {I + l)-encounter; only the link disappearing in the merger is drawn. Reference stretches in both structures thick 



Shrunken links: To get rid of the remaining contributions we consider an arbitrary structure involving a 2- 
encounter, and there focus specifically on the link ending at the first entrance port of the 2-encounter. Then where 
does that link start? The only possibility that we still have to consider is that it starts at the exit port of a different 
encounter (with an arbitrary number I of stretches). All other possibilities have been dealt with above: if the link 
would start at the exit port of the same stretch the corresponding original orbit would just include one stretch and 
one link; if the link would start at the exit of the second stretch, there would be a partner orbit with only one stretch 
and one link. Hence we must have a situation as depicted schematically in Fig. [S^ for I = 4. 
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Now let us consider a structure that looks topologically the same as Fig. [8^ but has the link just considered shrunk 
away. In that structure the two ports connected by the link in question (the black ports in the picture) are ignored, 
leaving altogether / + 1 entrance ports and I + 1 exit ports. These ports have to be connected as in Fig. [8K, with 
connections in the original orbits following the full lines and connections in the partner orbits following the dashed 
lines. If we leave out the black ports, this leads to connections as in Fig. [81d. We conclude that when the link between 
the 2-encounter and the ^-encounter is removed these two encounters merge into a single (/ + l)-encounter. The 
contribution of the new structure is proportional to the old one. Just the factor —-, -, r from the link and 

^ i(eA OP c+es or d) 

the corresponding factor i{eA or c + or d) from the removed entrance port are gone, and a change of sign is again 
incurred. Hence the two structures depicted have a chance to cancel. 

What remains to be done is a little bit of book-keeping (e.g., to account for the numbering of encounters, and 
for the possibility that removing links from different structures yields the same result). To do so, we single out one 
stretch in each structure as a reference; that stretch must not be the first stretch in its encounter. The number of 
eligible stretches (i.e. the number of all stretches L minus the number of encounters V) will be denoted by n. We 
now sum over all structures and over all possibilities of choosing references stretches inside a structure; to avoid an 
obvious overcounting we modify the summand of Eq. (j42p to include an additional divisor n. 

Now for each contribution where the reference stretch belongs to a 2-encounter we can find exactly one contribution 
with opposite sign where the reference stretch belongs to an i-encounter with I > 2. This is done by (i) merging 
the 2-encounter with a different encounter as described above, (ii) renumbering the encounters to make up for the 
loss of one encounter (i.e. decrease by 1 the labels of all encounters which are larger than the label of the removed 
encounter), and (iii) choosing the encounter stretch following the merged one (e.g., the fourth encounter stretch in 
Fig. [8)3) as our new reference stretch. The last step guarantees that the reference stretch ends up in any position but 
the first in its encounter, as demanded above. It is clear that the new contribution has the same n as the initial one 
(as both L and V have been reduced by one), hence the additional divisor 1/n weighing every contribution remains 
unchanged. 

Each contribution where the reference stretch belongs to an ^-encounter can be accessed by following this procedure. 
This becomes clear if we take a contribution as in Fig. [SJs and then undo the steps described above to arrive at Fig. 
[5^. The same contribution as in Fig. [SJs results from V different contributions as in Fig. [S^ since the removed 
2-encounter may have an arbitrary index between 1 and V . However this degeneracy is compensated by the change 
of the factorial in the denominator of Eq. (j42p . Since the remaining factors are the same but the sign is different, we 
thus see that every contribution where the reference stretch belongs to an ^-encounter exactly cancels with V "parent" 
contributions with the reference stretch in the 2-encounter. 

To summarize, we have seen that all off-diagonal contributions cancel due to two mechanisms: For quadruplets 
where at least one orbit is limited to just one encounter stretch (which is the first in its encounter) and one link, the 
structures where this orbit is assigned to different pseudo-orbits mutually cancel. For quadruplets where no orbit is of 
this type we singled out a reference stretch. We then distinguished between contributions where this reference stretch 
is part of a 2-encounter or an Z-encounter with I > 2. By the construction of Fig. [8] we were able to show that these 
two types of contributions mutually cancel. Hence for systems without time-reversal invariance we have = and 
only the diagonal term remains. 

VI. SEMICLASSICAL CONSTRUCTION OF A SIGMA MODEL 

An alternative proof for the cancellation of all off-diagonal contributions can be based on field-theoretical techniques. 
We shall show that the whole series (|42p can be summed to become an integral over a certain matrix manifold. That 
matrix integral turns out to coincide with the one known from the cr-model in RMT, and gives the same result as the 
diagonal approximation. 

A. Matrix elements for ports and contraction lines for links 

We recall that structures of pseudo-orbit quadruplets were defined by (i) fixing the numbers of encounters and 
stretches involved (the stretches being ordered according to the way they are connected in AU C and B U D), (ii) 
connecting the encounter ports by links, and (iii) distributing the resulting pre-reconnection orbits over the pseudo- 
orbits A and C, and the post-reconnection orbits between B and D. 

Now we encode each structure in (|42t by a sequence of alternating symbols Bkj and Bjk respectively standing for 
entrance and exit ports. We start with the first entrance port of the first encounter, continue with the first exit port 
of the first encounter, then list all further ports of that encounter, then of all further encounters. The indices of B, B 
indicate the affiliation to the pseudo-orbits. Namely, the second index of B and the first index of B determine whether 
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before reconnection the port belongs to A (case j = 1) or to C {j = 2). Similarly k shows whether the port belongs 
after reconnection to B (if k — 1) or to D (fc = 2). Due to the two possible values for both k and j we altogether 
confront four i?fcj's and four Bj^s. 

Our conventions for the numbering of ports imply that the ports corresponding to a B and the subsequent B 
(entrance and exit ports with the same number) are connected by an encounter stretch of an original orbit. Hence 
they must belong to the same original pseudo-orbit and have the same index j, as in Bk'jBjk- Similarly the ports 
associated to a B and the subsequent B (representing an exit port and the entrance port with the following number) 
must be connected by an encounter stretch of a partner orbit. Consequently they belong to the same partner pseudo- 
orbit and share the same index fc, as in Bj^B^j". The same applies for the last B and the first B in each encounter. 
In any case neighboring indices of successive port symbols must coincide. For a 2-encounter with two entrance and 
two exit ports as in Fig. [6] (a) we thus get the four-symbol sequence 

just like in a trace of a matrix product; note that the sub-indices of the indices j, k signal port or encounter stretch 
labels within the encounter. For a structure with arbitrarily many encounters indexed by cr = 1, 2, ... 1/ and stretches 
inside each encounter indexed by i = 1, 2, . . . l{a) the sequence of ports gives rise to the product 

V 

n (-S/C„l,i„l-Bj„l,fc„2-Sfc„2j„2 • ■ --Sj^n^j.^vi) ■ (45) 
cr=l 

Note that we do not imply summation over repeated indices either here or anywhere else in the paper. The indices j 
and k now carry two subindices each, cr = 1, 2, to label the encounter and i = 1, 2, . . . , l{a) to label the stretches 
within an encounter. Such "proliferation" of subindices notwithstanding there are still only four different -Bfcj's and 
four BjkS. Whenever possible below we shall suppress the "address labels" on k and j. If we take Bkj,Bjk to be 
elements of 2 x 2 matrices B = (^^l^^ bII^ ^^"^ similarly for B we may consider the expression above as one of the 

summands in the expansion of the trace product J1ct=i (-B-B)'''^^ 

The links of a structure will be indicated by drawing "contraction lines" between the respective exit and entrance 
ports. Each such line connects one entrance and one exit port of some orbit and its partner (recall that links do not 
change when shifting connections in encounters). The indices j, k of B and B indicating accommodation in original 
and partner pseudo-orbits therefore have to coincide. 

Example: In our above example of a product of four symbols the ports can be connected by links in two ways: 

I ^ ] I 1 I 1 

BkijiBj^k2Bk23iBjiki or BkijiBj^kiBkij2Bj2ki (46) 

note the agreement of indices of the connected ports. We can even complement these diagrams to topological pictures 
of the orbits. If we draw lower lines to indicate how the ports are connected inside the encounter we obtain 



I ^ [ I ^ [ 

Bki3iBjik2 fiajifi ifci and BktjiBj^k2 Bk2jiB j^kt 

or 

I J I 1 I 1 I 3 

BktjiB j^ki ■ Bfcij2-5 j2fci and ^fcui^jifci-Bfcijs-B jafei • (47) 

Here the upper pair represents the structures in the list (j43p and pertains to Fig. [SI The first picture shows the 
single orbit before reconnection (with encounter stretches connecting the ports as BB). The second picture depicts 

the two orbits after reconnection (with encounter stretches connecting the ports as BB). Analogously the lower pair 

represents the structures pertaining to Fig. [6](c). 

Diagrammatic rules: We can now translate the previous diagrammatic rules for the contribution of a structure in 
(I42p to the new language of contracted sequences of matrix elements: 

(i) The ports Bkj and Bj'k' connected by a link must have the same subscripts j = j' and k = k', and the 
contribution from that link in equals [—i(£AorC + or d)]~^ ■ We can thus write the factor provided by a 
contraction between Bkj,Bjik' as [— i(ej -I- £'k)]~^Skk'Sjj' where ej stands for (j — 1) or ec (j — 2); similarly e'f. 
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equals eB for fc = 1 and e/j for k = 2. (The primed energy will always refer to the post-reconnection pseudo-orbits B 
or D.) 

(ii) The indices in an encounter factor i {eAoiC + esoru) in (|42p are determined by the pseudo-orbits containing 
the first entrance port of the encounter. In (j45p the first entrance port of the ath encounter is denoted by 

such that the respective encounter factor can be rewritten as i(ej^j + ^fc i)- 

(iii) Finally Eq. ([55]) implies that there should be a factor (— l)"c+nD depending on the numbers of orbits included 
in the pseudo-orbits C and D. 



B. Wick's theorem and link summation 



We now want to sum over all structures, which in particular involves summation over all possible ways to draw 
links alias all contractions. A useful tool for that summation is provided by Wick's theorem. Consider the integral 



dbdb*e'""' f{b,b*) (48) 

over the complex plane with dbdb* = '^R^yi'"'' and /(6, 6*) a product involving an equal number of 5's and 6*'s, and 
Ime > to ensure convergence. Wick's theorem states (i) that such Gaussian integrals can be written as sums over 
"diagrams" where each b in f{b,b*) is connected to one b* by a "contraction line", and (ii) that for complex b,b* 
these diagrams can be evaluated by performing the Gaussian integral over contracted elements as if the rest of the 
integrand were absent. The contraction lines can then be eliminated as in 

J dbdb* bb*g{b,b*)e''''''' = dbdb* g{b,b*y^^^' ■ (49) 

Using this contraction rule one can stepwise remove contraction lines and the 6's and 6*'s connected by these lines 
and thereby in each step earn a factor — For complex b and b* Wick's theorem is just a fancy way of writing down 
the integral / dbdb*^^^^'' {bb*Y = pi^yrrFT; the factorial gives the number of possible ways to draw contraction lines. 

But remarkably. Wick's theorem also holds, up to a sign in Eq. ([49]) , if we replace b and b* by Grassmann 
variables rj and 77*. Grassmann variables are formally defined as anticommuting, rj-q* = —rfrj, and with integrals 
J dr] = J dr]* = 0, J drjrj = J dr]*ri* = 1 "^l. While 77 and rj* should in general be seen as independent variables, for 
the purposes of this paper it is convenient to imagine rj* to be the complex conjugate of rj and define the complex 
conjugate of 77* to be (77*)* = —7/. Imagining exp(ier7?7*) Taylor-expanded we obtain 

dT]dr]* 7^7* e'"'"' = J dT]dr]* 7777*6'"""' ^ = J ^vdv* e''"^'^' ■ (50) 

This is the analogue of (149]) with g = 1 which is the only case of interest since all powers of the Grassmann variables 
higher than 1 vanish. Below we shall apply contraction rules of the foregoing types to integrals over an arbitrary 
number of pairs of commuting and anti-commuting variables. Following traditions of the physical literature we shall 
refer to the anti-commuting variables as "Fermionic" as opposed to the commuting "Bosonic" ones. 

We now want to use Wick's theorem to represent sums over all possible ways of drawing contraction lines between 
Bfcj's and BjkS, for general products of the form This can be done if we declare these symbols to be either com- 
plex Bosonic or Fermionic variables and make the variables Bkj and Bjk between which we want to draw contraction 
lines mutually conjugate up to a sign, i.e., we set Bjk = iB^^. (The choice of the Bosonic or Fermionic character of 
the variables and of the sign will be made later.) Then all possible contractions are generated by taking the corre- 
sponding expression without contraction lines, and integrating with a Gaussian weight that has a term proportional 
to BkjB'^j in the exponent. The prcfactor of BkjB^- should be chosen as \{ej + e'^). Then Wick's theorem yields 
the desired factor [— i(ej -1- efc)]~^ for each link belonging to the pseudo-orbits j and k (up to a sign factor arising for 
Fermionic variables). If we also include the encounter factor i(ej^^ -I- e'j, ) we are led to expressions of the type 



± j d[B,B]e^^^.''=^.^'^''+''''^^>''^'^^ jj]^ ( 



±i(^J<.i + 4„i)-Bfc„ij,.i^j<,i,fc„2-Sfc„2j„2 ■ ■ - Bj^n^-jM^ijj (51) 
where the integration measure involves the product of all independent differentials, 

d\B, B] = X^d (RcSii) d (ImBii) d (ReB22) d (ImB22) dBi2dBUdB2idBX. . (52) 



All (convergent) expressions of this type contain all relevant structures with the appropriate link and encounter factors. 
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C. Signs 



The most delicate task is to fix the signs in ((5T|) and the Bosonic vs. Fermionic nature of the variables so as to 
produce the sign factor (— l)"c+nD ^j^g sum over structures (|^^ . We shall make the following choices which are 
inspired by the supersymmetric sigma model of RMT (see 0] and Appendix [B] ) but will here be justified purely 
on semiclassical grounds: We take -611,-622 Bosonic and -612,-821 Fermionic, and let the supermatrices B and B be 
related as 

where = (0-1) ^^'^ adjoint of a supermatrix is B'' = gV^^- Then we pick the signs in the Gaussian 

and the product of -B's and -B's of ((5T|) such that these terms can be written in terms of supertraces StrX = tr (t^X. 
We also add an overall minus sign. 

For the Gaussian integral we thus write 

- J d[B, B] exp (i(ei + e'i)BiiBi*i - i(ei + e'^)B2iB*^ + i(e2 + e'i)Bi2Bi*2 + i(e2 + £'2)52253*2) {■ • ■} 

= - j d[-B, -B] exp (^iStre'B-B + iStre-Bs) {...} = (({...})); (54) 

here the curly brackets {. . .} have the same content as in (jSip : to compact the notation we employ the matrices 

e = diag(ei,e2) , e' = diag(e'i, e^) ; (55) 

moreover, in the last member of the foregoing chain of equations we abbreviate the -B-integral with the Gaussian 
weight (including the negative sign factor) by double angular brackets ((. . .)). Since the Gaussian is not normalized, 
((...)) is not an average, i.e., ((1)) ^ 1. Instead we have 



(1)) = - j d[B, B] exp (i(ei + e'i)Bii-B*i - i(ei + e^).B2iB*i + i(e2 + e'i)Si2S*2 + i(e2 + £2)522^2*2) 

^ (ei +e2)(£2 +ei) ^ (e^ + e-D)(ec + eg) 
" (ei + e'i)(e2 + 4) ~ i^A + eB){ec + <^d) ' 



(56) 



reminiscent of the diagonal part ^d|ag '^^ Z^^\ The factor ((1)) remains after all contractions are performed and thus 
all B, B removed. Reassuringly it accounts, once multiplied with e'('''+'^^~'^c'~'^°)/^, for all diagonal quadruplets. 

When writing the product of B 's and B 's in (|5ip in terms of a supertrace, we still need to keep the indices fixed (note 
that no summation over port labels is implied!). We thus employ projection matrices Pi = diag(l, 0), P2 — diag(0, 1) 
and write (jSTj) as 

|n + eL) StrPfe„,SP,„,M„. . ■■Pj^u.^Bjl . (57) 

Note that just like for a trace, it is possible to cyclically permute the supermatrices appearing in this supertrace. 
If we also install the weight -^iy of a structure and the Weyl factor, we obtain 

„i(eA+eB— EC-e-D)/2 

zW(l/,{/(a),>„fc..}) = 

X I n + 'li ) Pk., BP,^, BPk^, . . . B 1 (58) 

which gives the cumulative contribution to Z^^"^ = Z^l^(l + zj^g) of all structures with V encounters, each having a 
fixed number of stretches /(a), a = 1 . . .V and a fixed distribution of ports {jcri, ka-i} among the pseudo-orbits; each 
allowed way of drawing contraction lines makes for one structure. The purely diagonal additive term is due to = 0. 

The only thing still to be shown for a proof of this statement is that our conventions incorporate the correct sign 
(— 1)""^+"" for each structure depending on the number of orbits included in C and D. To fill this gap we reexpress 
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the contraction rules (|49p and ([50)) in terms of B and B. Let us consider one way of drawing contraction lines in 
Eq. (|58p . and moreover select a particular contraction, either between ports within the same encounters (inside the 
same supertrace) or from different encounters. Then we can prove two identities connecting the contribution of a 
structure with that of another structure containing two ports and one link less (see for details of the derivation), 

Str {PkBP,Y)Str{XPyBPk') ■■}) = - .f'''^'"' ((Str (PkXP.Y) ...)), (59) 



Str [PkBPjUPyBPk'V) ■••]) = -^^^^ ((Str (^it^)Str {PkV) . . .)) . (60) 

In these formulae, contraction lines other than the selected one are not drawn and are assumed to coincide in the left 
and right hand sides. In terms of periodic orbits, the rules ([59l[60|l respectively account for the two ports connected by 
a link belonging to either different encounters or to the same one. The symbols X, Y, U, V represent matrix products 



as in dSll) with Y = B . . . B, X = B . . . B; onthe other hand, U ^ B ...B, V = B ...B. The terms S^j^Skk' i^iqiFT 
are consistent with the diagrammatic rules stated above: Links can only connect ports associated to the same original 
pseudo-orbit j = j' and the same partner pseudo-orbit k = k' , and contribute factors — j^^-qip-j . The contraction rules 
can be used to stepwise simplify our expressions for the structure contributions. Each step removes two matrices or, 
in other words, a link connecting two ports from a periodic-orbit structure, without changing other links. 

At first sight, the rules (|59I60[) do not seem to yield any sign factors. But we must pay special attention to steps 
where the contraction line pertaining to the last link of an orbit is removed, after removal of all other contraction lines 
of some "parent" orbit in previous steps. The B and B so connected represent the only ports of the remnant orbit. 
That orbit is still periodic, and therefore the two ports are connected not only by the link but also by an encounter 
stretch. The two matrices i3's and -B's then have just one projection matrix in between. If the orbit in question is an 
original one, the matrices follow each other like BPjB. We thus have to use contraction rule ([60]) with U = 1. Then 
the first supertrace on the r.h.s. turns into StrPj which equals 1 if j = 1 (the orbit belongs to pseudo-orbit A) and 
— 1 if j = 2 (the orbit belongs to C). As desired a sign factor —1 thus arises for each orbit included in C. 

An analogous result holds for partner orbits. Imagine that we remove the last contraction line associated to a 
partner orbit. Then the ports connected by this line must also be connected by a partner encounter stretch. The rule 
(1501) then applies with V ^ 1. The second supertrace on the r.h.s. thus equals StrP^ and yields — 1 if fc = 2, i.e., if 
the orbit belongs to D. 

To conclude, by turning B and B into supermatrices, we have successfully incorporated the sign factor (— l)"c+nD^ 
Other ways to deal with that sign factor (the replica trick) are discussed in Appendices [C] and |D] 



D. Towards a sigma model 

To determine Z^-^^^ we sum over all structures of pseudo-orbit quadruplets. In Eq. ([55]) we had already collected all 
ways of connecting ports with links. Now we sum over all possibilities of assigning ports to pseudo-orbits, i. e., over 
the indices j^i, k^i in (|58p . This means that we have to drop the projection matrices. To keep the factors i{tai + ejji) 
connected to the indices in the supertrace we use the diagonal offset matrices e, e' and write 



zW(T/,{Z(a)}) = 



i(eA-(-eB-cc-CD)/2 
V\ 



J]iStr U{BBy'^^^ +e'{BBf "'^ 



(61) 



Here the term Str e'(i3i?)'*^'^' contains the factor ej.^^ from ([58)) while Str e(i?i3)'('^^ includes tj^^. It remains to sum 
over the number V of encounters and over their sizes Z(cr), u = 1, 2, . . . , y. We let V run from to oo, the term V — Q 
taking care of the purely diagonal contribution, and obtain 



CX3 ^ 



^ v\ 

V=0 



CX3 

iStr ^ (i{BBy + e'{BB) 



1=2 



_ gi (eA+eB-ec-er))/2 



cxp 



iStr J2 UiBBy +i'iBBy 



1=2 




(62) 



Now recall that (j62p was derived using a semiclassical approximation, which converges only with large imaginary 
parts rj for all energy offsets. In the limit rj ^ 1 the _B-integral draws overwhelmingly dominant contributions from 
near B = whereas other B only make for exponentially small corrections of order . Therefore in the important 
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integration region we have ||-Bi3|| <C 1 and the geometric series in the foregoing exponent can be replaced by its sum. 
Upon explicitly writing the Gaussian integral we get the compact result 



(1) 

d[B, B] exp 



i (eA+CB-ec-eD)/2 



i Str e ■ 



BB 

1 - BB 



L 



d[B,B] exp 



i ^ ll + BB ^,1 + BB 
- Str I e ^ + e' 



1 - BB 



1-BB 




(63) 



The integration domain C must include the vicinity of the point B — B — and obey 1 — BB > 0, to exclude regions 
where the integrand is exponentially large; otherwise it is arbitrary. Under these limitations and 77 S> 1 the asymptotic 
1/e-expansion of the integral ([55)1 coinciding with the semiclassical series (|61l) is created by the saddle point of the 
exponent at B = B = 0. 

If we evaluate (|63|) in the limit 77 ^ 1 (see Appendix [B] ) we obtain the same result Z'^^-' as in the diagonal 
approximation, Eq Then the same 

Z(2) 

as in the diagonal approximation is obtained using the Riemann-Siegel 
relation Z^^\eA, es, ec, er>) = Z^^\eA, es, —£d, — ec)- We see once again that all off-diagonal terms in the generating 
function cancel. 



E. Comparison with the RMT sigma model 

In the so called rational parametrization of the zero-dimensional sigma model of RMT (see, e.g.. Appendix [B] ) the 
exact generating function Z is given by the integral 

Z{€A,eB,ec,£D) = - d[B,B]exp 
J Co 

identical in appearance with (j63p and with the integration measure and the structure of the supermatrices the same 
as there. However there are two important differences: First, the offset variables eA,B,G,D are all allowed to be 
real. Second, the complex Bosonic integration variables are assigned specific integration domains, the open unit disk 
l-Biil < 1 and the full complex plane I-B22I < 00. It is easy to show that one of the eigenvalues of BB is then between 
and 1 while the other one can take any negative value. Consequently 1 — BB > and the integration domain Cq is 
in fact a special case of C in ([63)) . 

Now suppose that instead of calculating the integral exactly we consider its high-energy asymptotics. We then get 
a sum of contributions of the two existing stationary points of the exponent. The first stationary point is at i? = 0; 
its contribution coincides with our semiclassical result Z^^K 

The second stationary "point" associated with the so called Andreev-Altshuler saddle corresponds to Bn = 0, B22 
00 (see Appendix |B] ). For 77 3> 1 its contribution is exponentially small compared to the first one and can be 
dropped, but it is of the same order for small or zero 77. Now importantly this contribution is related to Z*^^-* 
by Z^'^\e A, €3,^0,^0) = Z'-^^eA,eB,—^D,—^c) (see Appendix [B|) . This is the same relation as between the two 
summands in our semiclassical approximation of the generating function, based on the Riemann-Siegel lookalike 
formula and thus unitarity. As anticipated the generating function and the two-point correlation function of fully 
chaotic dynamics without time-reversal invariancc thus agree with those derived from the Gaussian Unitary Ensemble 
of RMT. 

The following extension of the comparison with the RMT sigma model requires some more knowledge of that model. 
Previously uninitiated readers may want to first consult with introductions provided in Appendix [B] and in Q or with 
the pioneering work of K. Efetov [4] . 

Q-integral and geometrical interpretation: We first compact the RMT matrix integral (j64p to 

Z = J dQe^'^'"^(° -0 ; (65) 

here the 4x4 supermatrix Q = iTA^^^T^^ is related to a diagonal matrix (called "standard saddle") A^^^ = 
diag(l, 1, — 1, — 1) by conjugation with the supermatrix T = ^J;^^- The matrices Q and A'-'^' obey the nonlin- 
ear constraints —Q^ — (A*^^))^ = 1. The integration domain for the Bosonic elements in Q can be shown to be the 
two-sphere S2 and (one of the disjoint halves of) the two dimensional hyperboloid H2. In particular, the complex 




BB ^,1 + BB 
1-BB 



(64) 
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variable B22 in the foregoing rational parametrization can be identified as a stereographic projection variable for the 
sphere ^2. The point B22 = corresponds to the "south pole" while the antipodal "north pole" is approached as 
B22 00. The complex variable Bu, on the other hand, is associated with H2- 

The foregoing rational parametrization of the manifold of matrices Q is just one of many. In fact, the manifold in 
question can be characterized without specifying any "coordinates" , and we would like to just state that character- 
ization here, referring to Appendix [Bl of the present paper for some details. For a geometrical interpretation of the 
matrix integral it is well to reorder the lines and columns in the 4x4 matrices Q,T,A as 1234 1324, so as to be 

able to employ 2x2 blocks as Q = (q^b Qff) ' Bose-Bose block Qbb and the Fermi- Fermi Block Qff contain 
only commuting entries while all anticommuting entries are assembled in the off-diagonal blocks Qbf,Qfb- In that 
"Bose-Fermi notation" (the previously used ordering is known as the advanced-retarded notation) 's'l the integration 
domain for the blocks Qbb, Qff in the sigma model integral become 

Qbb G i72 = U(1,1)/U(1) X U(l), Qff G ^2 = U(2)/U(l) x U(l), (66) 

where U(2) is the 2-dimensional unitary group, and U(l, 1) is the group of 2-dimcnsional pseudo-unitary matrices, 
i.e. the group of matrices obeying Wa^U = a^, where = diag(l,— 1). Geometrically, H2 and S2 define the 
two-dimensional hyperboloid and the two-sphere mentioned before. 

Andreev-Altshuler saddle point fs^l . 132 . Isdi l: We now comment on another specific parametrization of the RMT 
sigma model connected to the contributions Z'-^^ in our semiclassical approach. This parametrization is obtained 
from the above rational one by the transformation Q = RQ' R with R a permutation matrix interchanging rows 
and columns as 1234 —f 1432 and obeying R^ — 1. The pertinent Jacobian is unity and the integration manifold is 
unchanged. The exact generating function can thus be written as the "Q'-integral" 

Z = y"dQ'e^^*' '^''^(° -0^ (67) 

and the rational parametrization chosen for T in Q' = TA^^^T"^. We have thus recovered the Weyl symmetry, 
Z{eA,<^B,^c,^D) — Z{eA, ^B, —^D, —^c), since conjugation of the matrix of energy offsets with R just amounts to 
ec "CD- On the other hand, it is most illuminating to check that the "point" B = B = for Q' corresponds 
to Q = RA'-^^R = A(^\ the so called Andreev-Altshuler saddle. The two diagonal "saddles" are not within reach 
of one another through the rational parametrization with any finite B, B. Rather, A'^^ can only be approached as 
TA(i)r-i for B,B ^ 00 and vice versa. In particular, south and north pole of the sphere S2 are swapped in the 
Q'-integral relative to the Q-integral. In that sense one may say that the Q-integral (|65p rationally parametrizes Q 
w.r.t. the standard saddle while the Q'-integral ([S7|) rationally parametrizes Q w.r.t. the Andreev-Altshuler saddle. 
We would like to emphasize that either variant of the rational parametrization yields, upon exact evaluation of the 
superintegral, the full generating function; the inaccessibility of the antipodal saddle does not matter since all of its 
neighborhood is accounted for. 

As already stated, our semiclassical construction of the [B, i?]-integral corresponds to choosing the rational 
parametrization according to the Q-integral in since that choice for the RMT sigma model integral does yield Z^-^'' 
in a Gaussian approximation around B = B = 0. Had we done our semiclassical analysis for Z(^^(e^, eb, eci ££>), we 
would have ended up with a [B, i?]-integral corresponding, in the RMT context, to the alternative rational parametriza- 
tion in the Q'-integral in (p7|) : The Gaussian approximation about B = B = in ([57)1 does yield Z'^' (e^, cb, ec, ed), 
as is easily checked. In other words, the asymptotic contribution which stems from the infinitely remote integration 
point in the parametrization using the standard saddle, becomes the contribution of the stationary point B = B = 
in the Andreev-Altshuler parametrization and vice versa. 

VII. TIME-REVERSAL INVARIANT SYSTEMS 

We now propose to generalize the above to time-reversal invariant systems. For them, the possibility to revert orbits 
(and their pieces) in time makes room for more constructive interference. We shall see that the contributions of orbits 
differing in encounters no longer cancel but rather generate terms of arbitrarily high orders in 1/e. The summation 
of all relevant pseudo-orbit quadruplets can again be done with the help of a sigma-model type matrix integral. 

A. Diagonal approximation and ofF-diagonal corrections 

For time-reversal invariant systems the orbits of {A, C) can not only be repeated in {B, D) but may also be reversed 
in time and only then included in {B.D). Thus if we again formulate the diagonal approximation in terms of an 
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exponentiated double sum over orbits (|26p . the relevant contributions to the exponent involve pairs of orbits which 
are either identical or mutually time-reversed. This doubling of orbit pairs leads to an additional factor 2 in the 
exponent, and the resulting term in Eq. ([28)1 needs to be squared. 

Moreover, we now have to consider encounters where some of the stretches are almost mutually time-reversed rather 
than close in phase space. For instance, Sieber-Richter pairs as depicted in Fig. [2] contain one encounter with two 
almost time-reversed stretches. If we change connections inside this encounter, the two links outside the encounter 
still look almost the same in configuration space, but one of them must be traversed with opposite sense of motion, 
and this is possible only if the dynamics is time reversal invariant. 

Since time reversal changes the sides where a stretch enters or leaves an encounter, it becomes cumbersome to work 
with the notion of entrance and exit ports. Instead we arbitrarily refer to the ports on one side of each encounter as 
"left" ports and those on the opposite side as "right" ports. Encounter stretches may lead either from left to right or 
from right to left. The links may connect left and left, right and right, or left and right ports. 

Our formal definition of structures is then changed as follows: SI: When numbering ports we now make sure that 
in A U C the i-th left port is connected to the i-th right port, whereas in B U it is connected to the [i — l)-st right 
port. S2: When connecting ports through links we have to allow for arbitrary connections between ports on either 
side of the encounters. S3: Apart from dividing the orbits before reconnection between the pseudo-orbits A and C 
and the orbits after reconnection among B and D we have to choose the sense of traversal of every orbit. 

Again, several structures can describe the same quadruplet. In particular, we decide arbitrarily which side of each 
encounter is chosen as left or right. For V encounters there are 2^ equivalent possibilities of labelling the sides. To 
avoid overcounting, this factor needs to be divided out. 

The rest of our argument proceeds in analogy to Section IIVI Our semiclassical generating function turns into 

^ ^i(eA+eB-ec-eD)/2 f (^^ + £c)(£C + ^b) 

'^'"^ \{eA + eB)iec + €D) 

7(1) _ _J_( T\nc+nD Ilcnc i or C + £B or £)) 

^off - 2^ WV\^~ ' ff TTiv TZ Vi' ^ ' 

s^rcs^ V\ llii„ks(-l(eAorC+eBorD)) 

where the subscripts in the encounter factor depend on the first left port. 

We note that some of these structures already appear for systems without time-reversal invariance, or coincide 
with such structures up to the directions of motion in some of the orbits. One can show that the proof of section fVl 
generalizes accordingly, i.e., all structures trivially related to those in the unitary case mutually cancel. 

Let us use the semiclassical expansion (j68p to calculate the leading-order oscillatory contribution to the correlation 
function. The oscillatory term is obtained by application of /deAdes to 

(e^, eB,ec, eo) = ^^'^ (cA, es, -eo, -ec) (69) 



{eA - ec) (es - <^d) 



(eA + es) (ec + ed) 



-1 2 

(1 + Z^ff (eA, es, -en, -ec! 



and then taking the limit e^, es, ec, ^ e. The factor [(e^ — ec) (e^ — ^d)]^ tends to zero in this limit even after 
taking the derivatives. Therefore the oscillatory contribution to the correlator is purely off-diagonal. It must be 
connected to summands in Z^g (e^., es, — eo, — ec) which tend to infinity like [(e^ — ec) (es — e^)] ^, since multipli- 
cation with such terms and differentiation w.r.t. eA and eb is needed to eliminate the factor [(e^ — ec) (e_B — ^d)]^ ■ 

The respective contributions in Z^^ {ea, es, ec, £d) behave like [(e^ + £d) (ec + es)] . Since inverse powers of the 
energy offsets can only originate from the link factors, such terms can only be due to quadruplets that contain a link 
belonging to the pseudo-orbit A before and to D after reconnection, and another link belonging to C before and to B 
after the reconnection. In particular this means that the corresponding quadruplets must contain at least two orbits, 
both before and after reconnection. 

In the lowest order these are the quadruplets depicted in column 2:2 of Fig. [31 (We omitted quadruplets containing 
two copies of Fig. [5] or one copy of Fig. [5] and one Sieber-Richter pair. These quadruplets are irrelevant since the 
contribution of a disconnected diagram can be written as the product of its components, and the factor due to Fig. 
[5] vanishes by the argument in section FlVEn Among them the only diagram that explicitly requires time-reversal 
invariance is the one consisting of two Sieber-Richter pairs (the middle diagram in column 2:2 of Fig. [3]). In order to 
satisfy the above conditions on membership in pseudo-orbits we only need to consider the situation where one of these 
Sieber-Richter pairs involves one orbit included in A and one orbit included in D, whereas the other pair involves one 
orbit included in B and one orbit included in C. There are 2'^ — 16 ways to choose directions of motion in each of 
the four orbits involved, and two choices for considering either of the two encounters as the first, leaving altogether 
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32 equivalent structures. Their total contribution to {eA, es, ecj £n) is 

1 [i(eA + eD)][i(ec + eB)] 4 



32 X 



222! [_i(e^ + e^)]2[„i(,^ + g^)]2 



The leading oscillatory contribution to the correlator thus reads 



Cdouble-SR (e) 
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deAdes 



{eA - ec) (es - cd) 
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5(eA+CB+cc4-eD 



{eA + es)^ (ec + ^of 



2e4 



(70) 



which coincides with the predictions of RMT. 



B. Semiclassical construction of a sigma model 



To go beyond the leading order, we again write our semiclassical results in terms of a matrix integral. The derivation 
of this matrix integral is very similar to the unitary case. First we encode every structure by a sequence of symbols, 
each symbol standing for a port, and links depicted by contraction lines. Second, we replace every symbol by a 
Bosonic or Fermionic variable and integrate the product with a suitably chosen Gaussian weight. The contribution 
of every structure will then coincide with the corresponding contractions in the resulting Gaussian integral; the sum 
over all possible link connections will be provided by the integral itself. Finally, summing over all encounter sets we 
get an integral representation of the generating function known from the zero dimensional sigma model of RMT. 



1. Matrix elements for ports and contraction lines for links 



We again write down the sequence of ports starting with the first left port of the first encounter, continuing with 
the first right port, etc. The left ports will be denoted by B^k,p,j and the right ports by B^j^j^^k- The first and the 
second pair of indices of B respectively refer to the properties of the port before and after reconnection, and vice versa 
for B. As in the unitary case the Latin index j — 1,2 indicates whether the port belongs to the original pseudo-orbit 
A {j = I), or C (j = 2). The index k reveals whether after reconnection the port belongs to the partner pseudo-orbit 
B {k = 1) OT D {k = 2). 

The Greek indices fi, v account for the directions of motion through the port and the attached encounter stretches, 
both in the pertinent original and partner orbits: \i [i — \ the direction is from left to right in the original orbit; if 
/i = 2 it is from right to left. The index v =\,2 indicates the same for the partner orbits. The symbols B^k,fij, B^j^^k 
can be considered as elements of 4 x 4 matrices B, B; to simplify the notation we shall often substitute a capital Latin 
letter for the composite index like J = (/i,j). 

As for systems without time-reversal invariance the ports are ordered such that an element of B and the immediately 
following B represent ports that are connected by an encounter stretch of an original orbit. This means that they 
have to belong to the same original pseudo-orbit, and that the direction of motion (left to right or right to left) must 
be the same as well. Hence their indices J = {fJ-,j) have to coincide. Similarly each B and the immediately following 
B represent ports that are connected by a stretch of a partner orbit, and the corresponding subscripts K — {v, k) 
coincide. (When a B represents the last right port in its encounter, its subscripts v, k have to coincide with those 
of the first B.) Thus the arrangement of indices J and K is the same as in a product of matrices and we obtain a 
product of the same form as in Eq. (|45p but with j and k replaced by J, K. 

We now build in the links and depict them by "contraction lines" above the symbol sequences. Since for time- 
reversal invariant systems links can be drawn between ports on either side of the encounters, the contraction lines can 
now connect not only B to B but also B to B and B to B. Any ports connected by a link must belong to the same 
pseudo-orbit, before and after reconnection and hence their Latin indices must coincide, as in the unitary case. The 
Greek direction indices require new reasoning. A contraction between a B and a B stands for a link connecting a right 
port to a left port. In this case, if an orbit leaves one encounter at the right port it must enter the other encounter 
stretch at the left port. The direction of motion (left to right) is thus the same for both encounter stretches. The same 
applies for a direction of motion from right to left, and for original and partner orbits alike. Hence for contractions 
between B and B the Greek subscripts indicating the directions of motion must coincide, = fJ-2,t^i = 1^2- On the 

I 1 £ 1 

other hand, in the contractions B . . . BB . . . B the connected ports lie on the same side, and hence their directions 
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FIG. 9: In the presence of time reversal invariance a link can connect ports on the same side as well as on opposite sides of an 
encounter. The linked encounter stretches are traversed in antiparallel or parallel direction, correspondingly. 



of motion must be opposite. These possibilities are demonstrated in Fig. O We shall use a bar over i^, ^ to indicate 
that the port direction of the motion is flipped like 1^2 such that p, = 3 — fi. Again denoting the energy offsets as 
ei = eA, ^2 — ^1 — ^B, ^2 — have the link contribution 



I 1 

• [^i(ei+4)] Skk'Sjj>S^f,'Su,y' for B^k,fj.jBf^^j>yk' , (71) 

• [^K^i+efe)] Skk'Sjj'Sp,^'Si},y' for B^k,tJ,jB,j'k',f_i'j' and B^^^^kB^,' yk" ■ (72) 

The contribution of an encounter is defined by its first left port and can be written as i(ej„i + ^). And again, there 
will be a factor (— 1)"c+"d depending on the numbers of orbits included in C and D. 

2. From symbol sequence to Gaussian integral 

In order to represent additive terms in (|68p by contractions in a Gaussian integral we interpret the symbols 
Buk.pLj^Bfj^iji yk' as (Bosonic or Fermionic) variables. Again any two ports which may be connected by a link must 
be represented by a pair of variables which are mutually complex conjugate, up to a sign. Since rule ([7^ allows for 
contractions between two i3's differing in both their subscripts /i and v we thus need to have 

Byk,^J,j = ^B*,f.j^j . (73) 

To write the matrices B in a more elegant way, it is helpful to divide B into four 2x2 blocks, each block having 
fixed direction indices /i, and matrix elements indexed by k,j. The diagonally opposite blocks must be mutually 
complex conjugate up to a sign (which may even vary with j and k). We may write 



with the question mark on the equality a reminder of the as yet unspecified signs. The index p on bp, b* signals ^ = 
i.e., the ports associated with these subblocks are traversed in parallel directions by the original and the partner 
orbits. In contrast, the ports associated with 6o, 6* have ^ ^ v and their traversal directions by original and partner 
orbits are antiparallel. 

In addition, rule (|7ip allows contractions between matrix elements of B and B with identical subscripts. Hence 
these matrix elements, too, must be chosen as mutually complex conjugate (up to a sign). Since the subscripts in B 
and B are ordered differently, this means that B must coincide with the adjoint of B (again up to signs), i.e., 

B^B^ . (75) 

A Gaussian-weight integral of a product of mutually conjugate matrix elements will yield a sum over all permissible 
ways of drawing contraction lines. To obtain the correct contribution from each link the prefactor of B^k,vjB*^i^ in 
the exponent of the Gaussian again has to be chosen as i(ej + e'^,). Furthermore, due to the internal structure of B 
the sum over all B^^k,vjB*^^ ^- will include each pair of complex conjugate elements twice, which means that we have 
to divide the exponent by 2. Using the composite indices J,K and defining ej = tj and = independent of the 
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Greek part of J, K we can write the exponent in the Gaussian as J2jk ^'^i^J + ^'k)BkjB*j^j. Incorporating also the 
factor i(ej„i + j) for each encounter we obtain 

V 

X n (±i(ej<,i + ^'K^,)BK^^,j„^Bj^,jia2BK^„j„, ■ ■ ■ Bj^^^^^^k^,) ; (76) 

CT=1 

each (convergent) integral of this type includes the correct encounter and link factors for all structures with V 
encounters involving l{l),l{2) . . . 1{V) encounter stretches. 



3. Signs 

In order to incorporate the sign factor (— l)"c'+nr> g^j^j^ ^j^g correct prefactor of the diagonal approximation, we need 
to fully specify the matrices B and B: The off-diagonal entries in all four blocks of B are chosen Fermionic and the 
diagonal entries Bosonic, and the signs in (175)) are picked as 

here the multiplication with Uz in B means that each 2x2 block of B'^ is multiplied with az] taking the adjoint of 
our 4x4 supermatrices involves, apart from complex conjugation and transposition in the sense of ordinary matrices, 
a sign flip in the lower left (Fermi-Bose) entries of all 2 x 2 blocks. 

The signs in the Gaussian weight and the matrix product are again chosen such that both quantities can be written 
in terms of supertraces. The supertrace is now defined as the sum of the two upper left (Bose-Bose) entries minus 
the sum of the lower right (Fermi-Fermi) entries of the two diagonal blocks; in other words, the supertrace includes 
a negative sign for all diagonal elements associated to a Latin index 2 (regardless of the Greek index). The integral 
(j76p then acquires the form 

d[B,B] exp (^-Stre'BB + ^StreSsj (...) ee ((. . .)) , (78) 

with e,e' 4x4 diagonal matrices whose diagonal elements are ej = e^, ^'k — ^'k- After eliminating all contraction lines 
from any integral of this type we are left with the Gaussian integral 



(1)) = j d[B, B] exp ^^Str e'BB + ^Str eBB^ 

(ei + £',)^(£2 + e[)^ _ jsA + epfjec + £b)" 
" (ei + e'i)2(e2 + e'l)^ (e^ + eB)2(ec + enf ' 



(79) 



as in the diagonal factor for time-reversal invariant systems. Therefore, upon adding to (|78p the Weyl factor and the 
factor 1/{2^VI) from ([55]) we obtain the contribution to the generating function (incorporating the diagonal term) as 

„i(cA+eB-ec-«D)/2 

^ I n + JStr Pk^, BPj^, BPk^, . . . Pj^,,^,bJi , (80) 

where Pj [Pk) is a projection matrix whose diagonal element associated to J [K) equals 1 while all other elements 
vanish. 

We still have to show that with the choices made above the sign factor (— l)"c+nD £qj. ga,ch structure is obtained 
correctly. The contraction rules ([55]) and ([50]) we used to prove the corresponding result in the unitary case carry 
over. In addition it is now possible to draw contraction lines between two iJ's or two i?'s. These lines represent links 
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connecting two left ports or two right ports. The corresponding contraction rules are established in Appendix |3 they 
read 

Str {PKBPjY)Str [Pk'BPj.Z) . . .| = -^^—L-^S^k'^jj' ((str [PjYPkZ) . . .J) (81) 



Str {PKBPjYPK'BP.rZ) . . . j> = --^—L-^Skk'Sjj' ((str {PkYPjZ) ...)); (82) 



analogously for B. Here Y and Z are matrix products of the type required in ((80|) starting and ending with B. The 
sequences Y , Z are obtained from F, Z by (i) reverting the order of matrices; (ii) replacing all i?'s by B's and vice 
versa; (iii) replacing all projectors Pk by and Pj by Pj. Eqs. ([5T|) and (|82p again contain the familiar factor for 
each link. We just chose the more compact notation 

J=(Ai,j), J=(m,j), K^iiy,k), K={9,k) (83) 

The Kronecker deltas in the contraction rules (|81l82p imply that the connected ports must belong to the same 
pseudo-orbits and have opposite direction of motion. 

We can now use the contraction rules ([55]) . (|ST|) and (|82p and their counterparts for contractions of two i?'s 

to step by step remove the contraction lines corresponding to all links of a structure. In doing so we do not meet any 
sign factors apart from the steps where the final link belonging to an orbit is removed. This final link must always 
connect a left to a right port, since otherwise it would not be possible to close it to an orbit just by drawing an 
encounter stretch. Hence in this step we only have to use the rules ([5^ and ([50)1 for removing links between left ports 
(denoted by B) and right ports (denoted by B). By the same argument as in Section IVTl we can show that the final 
step yields a factor Str {Pj) for each orbit in AU C and a factor Str (Pk) foi' each orbit in B U D. According to our 
definition of the supertrace these factors are equal to —1 if the corresponding small indices are j — 2 (corresponding 
to the pseudo-orbit C) or fc = 2 (corresponding to the pseudo-orbit D). We thus obtain a negative sign factor for 
every orbit included in C or D, and these factors combine to the desired term (— l)"c+nD^ 



4- Towards a sigma model 



Summing over structures as in Section I VII and again using the assumed largeness of the imaginary part r] of all 
energy offsets, see we get the generating function 



d[B, B] exp 



/ LXJ \ . / LXJ 

^Streij2{BBy]+^Stre'iY,{BB) 



i 1 + ^^ ^/ i + BB 

- Str e ^ + e' 

4 \ l^BB l-BB 



^1=1 



(84) 



equal in appearance (up to the factor i) to (|63p for the unitary symmetry class. The matrix integral again sums up 
the contributions from all structures which form an asymptotic expansion in inverse powers of energy offsets. The 
integration domain in (|84p must include the vicinity of the saddle point B = B = 0, and be such that the eigenvalues 
of 1 — BB, 1 — BB be positive, otherwise it is arbitrary. 

To see that this asymptotic series agrees with the Gaussian Orthogonal Ensemble (GOE) of RMT, we can argue 
like in the unitary case: The sigma model for the GOE leads to a generating function as ([84|l but with a specified 
integration domain and without the restriction 3> 1. The large-e asymptotics of the integral is determined by two 
saddle points of the exponent both of which have to be taken into account unless 77 ^ 1. The so called standard 
saddle at i? = leads to a contribution dominated by small B and thus coinciding with ([84]) . The contribution of the 
Andreev-Altshuler saddle is related to Z'-^^ by Z'^' (e^, e^, e^, e^) — Z'^^^eat^b, —^d, —^c): again as in semiclassics. 
Our semiclassical results thus recover the asymptotics of the generating function in agreement with the sigma model. 

The same must then be true for C(e). The two-point correlation function of fully chaotic time-reversal invariant 
systems is thus asymptotically given by the sum of two divergent series in powers of e"^ in Eq. ([3]). 

Now we proceed to the full C(e). As is well known, an analytic function can, under rather general conditions, be 
uniquely restored from its asymptotic series by Borel summation [6|. That method involves term-by-term Fourier 
transform of the asymptotic expansion leading to a converging series, followed by the inverse Fourier transform of the 
resulting analytic function. The first stage of the Borel summation applied to the two components of the correlation 
function gives the spectral form factor K(t) both for r < 1 and t > 1. The inverse Fourier transform recovers the 
closed RMT expression for the correlation function in Eq. for all energies e, thus completing the semiclassical 
explanation of universality. 
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VIII. CONCLUSIONS 

We have shown that chaotic systems without any symmetries (beyond time reversal invariance) display universal 
spectral statistics. To obtain both the non-oscillatory and oscillatory contributions to the spectral correlator R{e) we 
had to (i) represent the correlator through derivatives of a generating function (ii) semiclassically expand Z as 
a sum over quadruplets of pseudo-orbits (sets of classical periodic orbits), and (iii) uncover a Weyl symmetry of Z. 
The non-oscillatory and the oscillatory part of the correlation function are due to two summands of Z related by the 
Weyl symmetry. The oscillatory one could only be recovered in a semiclassical approximation that explicitly preserves 
unitarity; in this approximation the oscillations are connected to a phase factor which accounts for the average level 
density. 

Constructively interfering contributions to the sum over pseudo-orbit quadruplets originate from orbits that are 
either identical or differ in close encounters in phase space. They define asymptotic ^-series. For dynamics without 
time reversal invariance those series terminate after the first term, both for the non-oscillatory and oscillatory part. 
For time-reversal invariant systems, the series are infinite ones but can be summed to give R{e) in closed form. 

To implement the summation we write the semiclassical series for the generating function Z as a matrix integral. 
The sum over all possible topologies then turns into an integral well known from the zero dimensional sigma model 
of random-matrix theory. 

For both the unitary and the orthogonal symmetry class, the correlator i?(e) comes out as the closed-form expression 
familiar from random-matrix theory. In contrast to the primary i-series which makes sense only for |e| 3> 1, the final 
result holds for any real value of the energy offset and yields, for e <C 1, the familiar power law i?(e) c>c |ep, with 
the exponent (3 revealing the degree of level repulsion characterizing each symmetry class. Our results thus imply 
the semiclassical explanation of universal level repulsion. Likewise, all other characteristics of universal spectral 
fluctuations based on the two-point correlator i?(e), like the so called spectral rigidity [s^, [13], are now semiclassically 
understood. 

The techniques here elaborated (in particular the analogy between orbit-based semiclassics and the sigma model) 
will be, and in fact are already helpful in tackling further challenges in the field. Examples are (i) localization effects 
in long thin wires [sH], (ii) non- universal effects related to the Ehrenfest time (which has no counterpart in RMT) 
[39| . (iii) a semiclassical treatment of arithmetic billiards [4ft] (where intrinsically quantum Hecke symmetries have 
strong action degeneracies as classical counterparts), (iv) spectral correlations of higher order [4l[ and the density 
of nearest- neig hbor spacings, (v) transport through ballistic chaotic conductors [H, \^ including the "full counting 
statistics" |43l | (vi) further symmetry classes, (vii) transitions between symmetry classes [43 ). and (viii) last but not 
least, the relative status of the ballistic sigma model and periodic-orbit theory [23j . 
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APPENDIX A: CONTRACTION RULES 

In this Appendix we want to derive the contraction rules ([50]) . (|5T|) . and ([5^ from Eqs. (|49p and ([50)1 . We 

start with systems without time-reversal invariance. First of all, contractions will be nonzero if the contracted 
matrix elements are mutually complex conjugate. Using the definition of the exponent according to Eq. ((54|) we 
obtain for contractions between Bf^j and (the latter agreeing up to the sign with S^/j/), 

Bi~B,,k'9{B, B)J^ = --^-J^S,,,Skk' {giB, B))) . (Al) 

Here Sj is a sign factor equal to 1 if j = 1 and to —1 if j = 2. To understand its origin one has to take into account 
the contraction rule for the Bosonic variables (j = k) yielding — and ([50]) for the Fcrmionic variables (j ^ k) 
yielding t^. An additional factor —1 emerges when j ^ k ~ 2 due to the negative sign in B22 — —-622 and when 
j = 1, fc 2 because of minus at the term proportional to B21B12 in the exponent of ([54]) . 
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Using Eq. (|Aip we obtain the rule (|59|) for contractions between B and B in different supertraces (inter-encounter 
contractions) which conies about as 



Str {PkBP,Y)Str {Pk'XP.yB) . . .J) = (j^tr {PjYPkB)Str {P,,BPk'X) 

I L 

SjYjkBkjSj'Bj'k'Xk'j' ■ 



= — —TT^jj'Skk' {{sjYjkSjXkj . . .)) 

= i ((Str (P.yP/cX) . . .» . (A2) 

Here we used the cychc invariance of the supertrace, then wrote out the supertrace in components, with the help of 
Str {PkZ) — SkZkk for any supermatrix Z, and finally used Eq. (|A1|) . 

Similarly Eq. (|60p for intra-encounter contractions between B and B follows from 



Str {PkBPjUPj.BPk'V) ...)) = (( SkBkjUj.yBj,k'Vk'k 



SkBkjBj/k'Ujj'Vk'k 



-jrSjj'Skk' {{skUjjVkk ■ ■ •)) 



= -T^^-rrSn'Skk' ((Str (P, (7)Str (PkV) ...)); (A3) 
H^i I ^k) 

here in the second line we could interchange Ujji and Bj'k', since a nonzero result arises only for j ~ j' in which case 
Ujj is Bosonic and thus commutes with all other variables. 

We now turn to time-reversal invariant systems. For these systems the derivation of rules ([59|) and ([60|) carries 
over directly, the only difference being that J, K are now double indices consisting of the Greek and Latin parts, 
J = (/i,j), K — (i^, fc); the sign factor sj has the values 1 for j — I and —1 for j — 2 regardless of the Greek 
index, and sk is defined analogously. In addition the fact that the matrix elements of B and B are pairwise mutually 
conjugate allows to draw contraction lines between two B's or two i?'s. We thus obtain further contraction rules ([5T|) 
and ([82]) for contractions between two B's, and analogous rules for contractions between two B's. To derive these 
contraction rules we use the important symmetry of B as defined in (|77p , 

= ~B Ei?*!]* = -B (A4) 



where B* is the transpose of B, and we have introduced Y. — ( i'^q) and = S ^ = i^^ qJ ■ (Readers who want to 

check (|A4p are reminded that transposition of a supermatrix involves interchanging the indices of the matrix elements 
and afterwards flipping the sign of the lower left elements in each 2x2 block.) Similarly conjugating the projection 
matrices Pr-, Pj we get 

EFa'S* = Pr , SPjS* = Pj ; (A5) 



Now let us consider the matrix products Y and Z on the l.h.s. of l|8ip and ([821) . involving alternating sequences of 
-B's and B's and interspersed projection matrices. Under transposition and conjugation with E these products have 
(i) the order of matrices inverted (due to the transposition), (ii) B replaced by B and vice versa (due to (jA4[) ). and 
(iii) the subscripts J and K of all projection matrices replaced by J and K (due to (|A5[) ). Since the numbers of i?'s 
and i?'s are equal the sign factors from (|A4p all mutually compensate. 

Thus equipped we can prove rule (j8ip for inter-encounter contractions between two B's. By using the invariance of 
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the supertrace both under transposition of its argument and under conjugation with E we write 
Str {PkBPjY)SU {Pk'BPj,Z) . . \ = (sty {PKBPjY)Sti [T.Z^Pj.B^Pk'T.^) . 



Str {PKBPjY)StT {ZPj,BPk. ) • 

! , / jV'^j?,/^' ((str [PjYPkZ] ...)). (A6) 



Here in the third Hne Z denotes the product obtained from Z by steps (i)-(iii) above. We have thus proven Eq. ([81 
We now move on to rule (1821) for intra-encounter contractions between two i?'s, as in 



SiY{PKBPjYPK'BP.j,Z)...Ji . (A7) 

First of all note that the contracted matrix elements are Bkj and Bk'J', the latter coinciding up to a sign with B*^ j 
(cf. Eq. (|77|)'). Hence the only relevant case is that J' — J, K' — K. For this case we now write 



// r 



Str [PkBPjYPkBPjZ] ...)) = ((^skBkj{YPkB)jjZjj, ...)) (A8) 



and evaluate the bracket in (IA8D as 



{YPj,B)jj = [E*(Si?*S*)(I]P^E*)(Sy*E*)E]*^^- 
= {ll^BPKYY^j _^ 

= sjBj^kYk.j. (A9) 

Here we first replaced the matrix sequence by itself double transposed; note that double transposition of a supermatrix 
(not an identity operation!) leaves its Bosonic elements like the one with the subscripts J J unchanged. Then we 
invoked Eqs. (IA4|A5|) and noted that for the matrix elements at hand transposition just amounts to interchanging 
the two subscripts. In the third line we used that only those elements of S, E* are nonzero for which one subscript is 
equal to another one barred, i.e., the port directions of motion are opposite. Finally, we used Ej jEj j = sj Inserting 
(|X9|) into (|M1) we now obtain 



Str {PkBPjYPk'BPj,Z) ••■])- ^skBkjsjBj^kYkjZj^j, . 



= \ , . hj'^K,K' ((str {PkYP-jZ) . . .)) (AlO) 

which is ((82|l . Here we used (jAl| . defined Y in analogy to Z, and in the end rearranged the result as a supertrace. 
One easily shows that rules analogous to (|5T|) and ([5^ also hold for contractions between two i3's. 



APPENDIX B: SUPERSYMMETRIC SIGMA MODEL 

We propose to briefly review the basics of the sigma model; details can be found in d, 0, [H, HH, [H, H^l • The sigma 
model is a way to implement averages over random matrices. The appropriate random-matrix ensemble for systems 
without time-reversal invariance is the Gaussian Unitary Ensemble (GUE), consisting of complex Hermitian N x N 
matrices H with a Gaussian weight proportional to e"'^*'^^ ; the matrix dimension N is sent to infinity. The pertinent 
generating function Z looks like (O, but with the energy average (. . .) replaced by the GUE average. To perform that 
average, three variants of the sigma model are available, all based on representing determinants by suitable Gaussian 
integrals. The supersymmetric variant deals with superintegrals over both ordinary and Grassmannian variables while 
two replica variants exclusively use either Grassmann or ordinary Gaussian integrals. 
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1. From GUE to supermatrix integral over saddle-point manifold 

In the supersymmetric approach each of the determinants in the generating function ([5]) is represented by a Gaussian 
integral, Fermionic for the numerator and Bosonic for the denominator, 

Z(eA,es,ec,er>) (^1)^ ^ II J ^ ^^"^ ^ ^^"^ d7]*(j^dr]cmd-qhmdilDn}j (Bl) 

-iE,„,„ {-z'A^[E+<iA-H)r^r.ZAr.+ZBr^{E-<iB-H)rr,r^ZB■n} 



at the expense of the overall factor (—1)^ the signs of the "Bosonic" quadratic forms z\{. . .)za and . .)zb have been 
chosen so as to make the Bosonic Gaussian integrals convergent, in consistency with lin.tA.B,c,D > 0. rjc, 11^,710,110 
are anticommuting Grassmann variables. As usual, the asterisk on commuting variables means complex conjugation. 
On anticommuting variables the asterisk may, for the purposes of the present paper, be given the same interpretation 
(complemented by (?/*)* = —rj, see [sj). 

At this point hooking onto the presentation in Chapt. 10 of Ref. 0] ^, we introduce the supervector 



^= II - 'f. (B2) 
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VD 
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\ lie 



where each of the four $q's has N components ^am] the ordering of the in (|B2p is known as the advanced- 
retarded notation. To accommodate the relative sign of the Bosonic quadratic forms mentioned and the different 
energy arguments above we employ the (4 x 4) matrices. 



L = diag(l, 1,-1,1) = 



1 



E = E^^_ diag(.., -ec) - - ^ (B3) 

whose block form will be convenient in the future. Our generating function may then be succinctly written as 

Z = (-l)^y" d^*d^ {^exp|i$ti[i/ ^ ^ (B4) 



with (i<i>*d$ denoting the integration measure written explicitly in (jBljl . The Gaussian average over 
Hamiltonian matrices H with a normalized Gaussian weight oc e~"2"*'^ leads to Z = (— l)^exp( — 
27v Sa,/ EaLfi J2m^n ^*am^an'^*0n^ I3ni) ■ The Hcrmitian 4x4 supermatrix 

N 

i 



TV 

allows to further compact Z as 



1 ^ 

^ ^ ^c.n.<^*f3m (B5) 



Z=(-l)^y exp|ArStr[-i^QL- i(QL)2]| (B6) 

where the symbol Str denotes the supertrace. The integrand would be invariant under the transformation Q — > TQT^ 
with pseudounitary matrices T defined by 

T^LT = L , (B7) 



We follow the notation of ^ apart from reordering the elements of $ and setting the parameter A in the random-matrix integrals equal 
to 1. 
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if the (real part of) the energy matrix were proportional to the unit matrix. Due to 

p{E = 0) = N/tt (B8) 

the energy offsets ca.b,c,d /'^•t^P are of only order which means that in any case this symmetry persists to leading 
order in 

The exponent in (jB6[) is quartic in the integration variables <i>, $*. Upon employing a suitable Hubbard-Stratonovich 
transformation Q we render the ^-integral Gaussian and end up with an integral over 4x4 supermatrices Q, 



Q'^ + 2\n{E~Q) \, (B9) 



r N 

dQ exp < — 2^^^ 

with a flat integration measure dQ. For a details of the transformation, in particular the overall minus sign and the 
integration range in (jB9|) . we refer the reader to Ref. 3]. Important for what follows is that the integration range is 
restricted to matrices Q diagonalizable by pseudounitary transformations T. 

The large- iV limit allows for a saddle-point approximation. The saddle manifold is determined by 

Q = ^J— . (BIO) 

E-Q 

We may confine the discussion of this equation to £' = 0, which means that we evaluate Z in the middle of Wigner's 
semicircle. As seen the energy offsets E — E are of order 1/N and thus irrelevant for saddles, we may even drop E 
entirely. The saddle-point equation then simplifies to = — 1. 

Next, we expand the exponent in Eq. (|B9p about the saddle manifold = — 1, to first order in the energy offsets 

and still sticking to the center of the Wigner's semicircle, E — 0, such that E = —jk^ ( o ■ ^i^*^^ the supertrace of 

the unit matrix vanishes we get Str [Q2 + 2 ln(£^-g)] = StT[\n{E^ -2EQ + Q^)] ^ StT[\n{-l) + \n{l + 2EQ)] « Str2£;Q, 
where we again neglected corrections of higher order in 1/N. Going to the limit N oo we get the integral 

dge^^'"^(» -0 . (Bll) 

=-1 

The eigenvalues of the matrices Q must be equal to ±i. This means that the saddle-point manifold can be divided 
into 5 coniugacy classes differing by the numbers n+, n_ of the eigenvalues -l-i and — i. In fact, only one of them is of 
relevance [3j, with n+ = n_ = 2. Since the matrices Q must be diagonalizable by pseudo- unitary transformations the 
relevant Q's can thus be represented as 

Q = iTA(^'T"^ (B12) 

with the diagonal matrix 

A^i) =diag(l, 1,-1,-1). (B13) 



2. Rational parametrization 

There is an ambiguity in the choice of the matrices T — ^""^ ^ . Due to the degeneracy of the eigenvalues of 
A^^\ Q remains unchanged if we replace T by the product TK with K any invertible block diagonal matrix. Choosing 
K = ^"^Ij" t"^) have 

y ' ^ ' y-Bil-BB)-^ (l-BB)-i J ^ ' 

with B = TarT-^^ and B = TraT-^ ■ The pseudounitarity of T entails the matrix {TK^LiTK) = K^LK = 
^{Taa )HTaa ) ^^.^-i^^^ (T^^^ block dlagonal. On the other hand 
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such that block diagonaHty yields B ~ cTzB'^ ■ With the latter relation between B and B established we can write 
these 2x2 supermatrices as 



B 



Bii Bi2 
B21 B22 



St 



Bii B21 

^12 -°22 



B = a.B'' = 



Bii B21 
B12* -B. 



22 



(B16) 



note that the diagonal elements -Bii,i?22 are Bosonic (commuting) variables whereas the ofF-diagonals Bi2,B2i are 
Fermionic, in accordance with the advanced-retarded representation we are using. From this point on we rename TK 
as T. 

The saddle-point manifold results as 



l + BB 



-2B(1 - BB)-^ 




l-BB 



The final integral over the manifold becomes an integral over the 2x2 matrices B,B, 
Z / d[B, B] exp 



AstrTA^T-if'' ° 
2 V -e 



- / d[B,B]exp 



2 \ l-BB 



e2 



(eA-l-«B-ec-eD)/2 



d[B, B] exp 




(B17) 



(B18) 



i Str e 



1 - BB 



BB 
1 - BB 



(B19) 



with the flat measure d[B,B] — ^-§^^^-^^dBi2 * dBi2dB2idB2i- The integration ranges for the Bosonic integrals 
over Bii and ,622 are respectively < 1 and |i?22p < 00. Note that the foregoing integral representation of Z'^^ 

equals the one established semiclassically, see 

We now proceed to the asymptotics of the remaining integral (jBlSp in the limit of large energy offsets, |eA,B,c,r>| ^ 1- 
To that end we employ one more saddle-point approximation®. The only finite stationary point of the exponent 
in (|B18|) occurs a.t B = B = 0. Let us assume, momentarily, that the energy offsets are provided with a large 
positive imaginary part, 77 3> 1. Then this will be the only critical point of the integrand responsible for the high 
energy asymptotics of the integral (jBlSp . To get an asymptotic series of Z in powers of we expand the matrix 
(1 — BB)~^ and (1 — BB)~^ in the exponent in powers of _B, B about the unit matrix and further expand the resulting 
exponential about its Gaussian part. Scaling the integration variables as B — B' j ^J\t\ and similarly for B with 
a typical dimensionless energy offset e close to tA.B,c,D we see that the higher-order terms in the expansion indeed 
become small perturbations in the high-energy limit. The Gaussian superintegral giving the leading-order contribution 
reads 



^(1) ^ _ gi(eA-l-ef3-ec-«r>)/2 
J2 0. j2 



(B20) 



d Bii d B22 ^UBr,Bii(eA+eB)+B22B:;,(ec+eD)} 
J dB^2^i2dB^idB2ie'^'^i2Sl2(eB-|-ec)+B2*lS2l(eA+eo)} 



The integration domain for the Bosonic variables can be extended to the full complex plane; in view of 77 3> 1 the 
error thus introduced is exponentially small. Upon doing the remaining elementary integrals we find 



Z{eA, iB,ic, ~ Z'-^\eA,iB, ec, e 



ec)/2 



{<^A + <^d){<^b + eg) 
{CA + e-B){f-C + (d) '' 



ry » 1. 



(B21) 



The previous saddle-point approximation had the matrix dimension A'^ as its large parameter and we could take tA,B C d/^ — > as 
A*^ — ► 00. No conflict arises with subsequently considering the large-e asymptotics. 
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3. Evaluation of the matrix integral 

In fact the result (|B21[) fully exhausts Z^^^: In agreement with the Duistermaat-Heckmann theorem [i^, all further 
terms of the perturbation series vanish for the unitary symmetry class; in the framework of the semiclassical theory 
this was already demonstrated in Section |Vl 

To prove the absence of perturbative corrections in the present RMT context we transform integration variables 
such that the integral in (|B18[) becomes Gaussian and similar to (|B20p . To that end we introduce an analogue of the 
singular- value decomposition of the matrix B, 

B = UBV'\ B = diag{xB,XF) (B22) 
with unitary U and pseudo-unitary V, and parametrize the latter matrices as 

V rj* 1-7777*72/ V &-">"=■) \ -r* l + rr*/2y ^ ' 

The matrix B_ is diagonal with positive (numerical parts of the) diagonal entries, with the Bose-Bose entry obeying 
< < 1, due to |i3ii| < 1. Obviously, B = VdzBU^^ = VBU^^ — diag (xb, — xf), with the Fermi-Fermi entry 
negative. It further follows that the matrix products BB and BB are respectively diagonalized by U and V, 

U-^BBU = V~^BBV = BB= (""^^ ^ 2) = diag (h, h), (B24) 

with the eigenvalues 0<ZB<l,Zi?<0. The matrices C/, V also diagonalize the rational functions of BB and BB 
appearing in the sigma-model integral 

The representation (jB22p can be inverted to express the variables on the right-hand side in terms of the original 
elements of B, B. The uniqueness of the decomposition of our B is thus ascertained. Moreover, the eight new 
variables Ib, If, 4>b, 0f, '7*, '''*, t can be employed as integration variables instead of the original elements of -B, B. 
A somewhat tedious calculation of the Berezinian yields the measure 

7rn ni dlBdlFd(j)Bd(t)Fdri*dridT*dT ,-r,^n\ 
d[B, B] = 4,.(;^„;^). • (B26) 

Our next step is yet another change of integration variables which will be elements of matrices analogous to B. B 
but with Ib^f replaced by mB,F, 

C = [7diag(V7^, 7=^)^-^ , C = a^C^ . (B27) 

We obviously have 

CC=^^, CC=^^, (B28) 
1~ BB l-BB 

such that the integrand in the sigma-model integral becomes Gaussian in terms of the matrices C and C. Moreover, 
in complete analogy with (|B26P the integration measure can be written as 

jr^ /=.! dmBdmFd(l)Bd(f>Fdrj*dr]dT*dT 

d[C,OJ = "2 . (,B29j 

47r^(7nB — mF) 

Comparing the measures (|B26l) and (|B29l) we get 

d[C,C]-d[B,B]ji^^^^/-^^^d[B,i3]. (B30) 
[niB — mF) diB diF 

Therefore, in terms of the integration variables Cik, C*^ the sigma-model integral becomes Gaussian. Again assuming 
that the energy offsets are all provided with the large positive imaginary part 77 we can extend the integration domain 
for the Bosonic variables to the full complex plane, at a risk only of corrections of the order exp(— 77). Then the resulting 
C-integral will become identical with the leading-order integral ()B20p . We have thus proven that contribution of all 
higher-order terms in the expansion of the integrand about the saddle point _B = cancels, and the asymptotics of 
the generating function, in all orders of is exhausted (for 77 ^ 1) by (|B2ip . 
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4. Andreev-Altshuler saddle point 

If the imaginary part rj of the energy arguments is small or zero Z'^-* does not capture the full high-energy asymp- 
totics of the generating function: an additive component Z^^-' appears due to another, the Andreev-Altshuler or 
"anomalous" stationary point B^^. To show it let us use the parametrization (|B22p ) introduced in the previous 
subsection; note that the elements of the diagonal matrix B satisfy < xb < I, < xp < oo . Using (jB24p . (jB28p we 
can rewrite the supertrace in the exponent of the sigma-model integral like 

Str (i + i' l±^] = {y-^eV + U-'e'U)^^ i±4 - {V-'eV + U-'i'U)^^ . (B31) 

\ 1- BB 1- BB ) ^ ' 1 - xl ^ l + xj. 

At a stationary point the derivatives of this expression by xb and xp must be zero, from which it follows that xb must 
be zero, whereas xp can be zero or infinity. Consequently there exist two stationary points in the matrix space, (a) 
the standard saddle B = giving birth to Z^^\ and (b) the so called Andreev-Altshuler saddle i?^^ = diag (0,00(7) 
where oo^ stands for the infinitely remote point of the complex plane. The anomalous saddle is responsible for the 
part Z'^^ of the generating function exponentially small when 77 ^ 1 but of the same order as Z^^^ if the energy offsets 
are real. 

An elegant way to find the contribution of the anomalous stationary point involves retracing the steps leading to 
the rational parametrization and finding the point Q^^ of the manifold (|B12p associated with B^^. Let us represent 
B^^ and B^^ as limits B^^ = lima^o PK~^, B^^ = hma^o PKZl with the diagonal 2x2 matrices 



for simplicity we take a real, i.e., we approach the infinite point of the complex plane along the real axis. When a is 
finite the respective 4x4 matrix T — f ) factorizes. 



T. ^ ( £ ) ( ^.-'^ ^\ ] . (B33) 



Since the second factor in Ta commutes with A*-^^ we find 

= i lim r„A(i)T-i = ii?A(i)i?-\ 



10 



_ f Ka P \ _ I 1 
10 
10 



R = R-={ P = 0010 I- (B34) 



The matrix R describes a permutation transposing the second and fourth entries, i. e., the Fermi-Fermi ones. We 
thus have 

A(2) = RA^^^R-^ = diag(l, -1, -1, 1) . (B35) 

To get the part of the high-energy asymptotics connected with the Andreev-Altshuler stationary point we repa- 
rameterizc the integration variables in (jBlip as Q' — RQR. The Jacobian of the transformation Q ^ Q' is unity. 
Denoting by E' the offset matrix with the Fermi entries transposed, E' = RER = diag(eB, — ec, — e^i, eu) we bring 
the Q-integral to an equivalent form Z — — Jqi2^_i dQ' e^^*'"' ^ . Now we can set Q' = TK^^^T~^ and rationally 
parameterize T which leads to a superintegral differing from (jBlSp by the interchange ep, <-*■ —ec- Expanding the 
exponent around the stationary point B = B = corresponding to Q' — iA^^\ 01 Q = Q^^ = iA^^^ we obtain an 
asymptotic series in which truncates on its first term, Z'^^^ — Z'^^^jg^^-cp. Adding the contributions of both 
stationary points we get the full GUE generating function Z — Z*-^^ + Z*^^'. 



APPENDIX C: BOSONIC REPLICAS 



The main technical difficulty in establishing a one-to-one correspondence between semiclassics and the sigmal 
model of RMT are the spectral determinants in the numerator of the generating functions. These determinants lead 
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to unwelcome sign factors (— l)"c'+nD jjj |;jjg semiclassical treatment, and in RMT they cannot be expressed through 
Gaussian integrals over complex numbers. In the supersymmetric treatment the way around these difficulties is to 
introduce anticommuting Fermionic variables. These variables incorporate the appropriate signs in semiclassics, and 
allow for an integral representation of non-inverted determinants in RMT. In the present appendix we want to present 
an alternative solution based on the so-called replica trick. (In Appendix [D] wc will consider a combination of the 
replica trick with Fermionic variables.) 

While the replica trick offers considerable technical convenience, there is a price to pay: First, the replica trick is 
hard to justify mathematically; second, at least in the present (Bosonic) variant only the part Z^^^ is accessible. That 
price has let replicas fall into some disrepute. From a semiclassical point of view, the price looks more affordable. 
We have seen Z^^^ to yield the full two-point correlator with the help of the "crosswise" procedure of Sect. Ill Al or. 
equivalently, Z*^^) to be accessible from Z^^") by the Riemann-Siegel lookalike of Sect. Ill Cl 



1. Semiclassical construction of a sigma model 

To construct a replica sigma model from semiclassics, we start from formula ([57)1 . We the assume that the pseudo- 
orbits C and D are each divided into r — 1 subsets yielding altogether 2r pseudo-orbits. We denote the pseudo-orbit 
Ahy j = 1, and the r — 1 parts of C by j = 2, . . . , r. Similarly B is indicated by fc = 1, whereas the r — 1 parts of D 
are labelled by fc = 2, . . . , r. Each of the nc orbits in C can be included in any of the r — 1 subsets of C, and similarly 
for D. This leaves altogether (r — ly-c+no -^^rays for distributing the orbits of C and D among subsets. If we now 
formally take the limit r 0, this number of possibilities turns into the factor ( — needed in our semiclassical 
treatment. We can thus drop the sign factor if we write zj^^ as 



7(1) M - V ^ ncnc»(£j+£fc) .(^.^ 



where we sum over structures of 2r-tuples rather than quadruplets of pseudo-orbits. This means that now every way 
of distributing the orbits in C and D among the subsets is considered to give rise to a separate structure. The indices 
of ej, e'f. in each link factor are chosen according to the pseudo-orbits j, k that the link belongs to; the corresponding 
indices in the encounter factor are chosen according to the pseudo-orbits the first entrance port belongs to. In line 
with the above definitions of j, k we take 

ei = eA, e2 ~ ■ ■ ■ = = ec, 

e'l-es, e^ = ... = e; = ei5. (C2) 

For the diagonal approximation we write 
Z\'^^ = limZ^iVr) 

diag j._^Q off V / 

^ e^E,^,/2+.E.eL/2 -Q (_^(e,.+4))-i (C3) 

j,k—l...r 

where the factors —i are just inserted for later convenience; they cancel mutually in the replica limit r — > 0. Now can 
now construct a sigma model similarly as in the main part: Each entrance port belonging to an original pseudo-orbit 
j (either A or a part of C) and a partner pseudo-orbit k (cither B or a part of D) is indicated by a symbol B^j. 
Correspondingly each exit port is denoted by a symbol Bjk- The indices of i?'s and B^s belonging to the same 
encounter have to coincide like the subscripts in the trace of a matrix product. We can again sum over all ways of 
connecting these ports by links with the help of a Gaussian integral over B^j, Bjk- This is done in the same way as in 
the supersymmetric approach, with only one difference: Since the sign factor (— l)"c+nD j-^^^g ^jeg^ eliminated by the 
replica trick, all sign factors in the previous treatment can be omitted including those in the commutation relation of 
B,B. The coefficients Bjk,Bkj can thus be taken as ordinary complex variables with Bjk = B^j (i.e. B = i?^). All 
supermatrices are thus replaced by ordinary matrices which are now larger due to the r different values of j and k. 

To go through the argument in detail, we split the sum over structures into sums over the possible numbers and 
sizes of encounters, assignments of ports to pseudo-orbits, and ways of connecting the ports through links: 

E - E E E ■ (C4) 

structures ^ cnc ports links 
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The three subsums are done step by step. 

Link sum: We start by summing over hnk connections for fixed encounters and assignment of ports to pseudo- 
orbits. We denote by Lj^, L^"' the numbers of entrance and exit ports forming part of the original pseudo-orbits 
j = 1, 2, . . . r and the partner pseudo-orbits k = 1, 2, ... r. The hnks now have to connect each of the -^^j?"* exit ports 
of j, k to one of the L^^ corresponding entrance ports, and we can conclude L^" = L^j^ for all j and k. With that 
condition satisfied, all connections between entrance and exit ports associated with the same j and k are permissible, 
i.e., there are Yijk ^kj^- ^ Yijk -^fe/*' possibilities. For each j, k we obtain a combinatorial factor 5^ in ^outL)."!. 

We now combine the latter combinatorial factor with the Li-i + 1 contributions ., , , . , one from each link within j 
and k and one from the diagonal approximation. We can represent the product as an integral over complex parameters 
Bkj, where we let Bjk — B'^j, 

(-1(6, +4))^(^^5+'^<5i.5,L-i^"! = |d^Sfe,e'fe+^^)l^^^l^(i?,,)^S(B,,)^S" (C5) 

with cPBkj — ^ fQi- each pair j and k. The integral converges due to Imej,Ime'j, > 0. 

The only terms not yet taken into account are the encounter contributions i(ej -f- e'j,), and the factor 
gidZj 'j+IZfc efe)/2_i^ from Eqs. (|C1|) and (jCSp . With these terms included the link sum yields the generating function 

Z''\r) = E E ^ fllK^. + 4)) n / d^S.,e'fe+^^)(Hl^..l^)(s,^.)^S(S*^.)"-' • (C6) 

#cnc ports ' \cnc / jk 

We have thus re-expressed zi^'^ as an integral over arbitrary complex r x r matrices B. 

Port sum: Next, we sum over all possibilities to assign entrance and exit ports of encounters to pseudo-orbits j 
and k. These assignments will determine L^" and L^^^ and thus the numbers of factors Bjk for the summands in 

Eq. (jC6|) . To perform the sum in a compact manner, we arrange the factors Bjk such that all factors originating 
from ports of the same encounter are grouped together. Each of these groups starts with the factor associated to 
the first entrance port, followed by the first exit port, the second entrance port, etc., up to the last exit port. This 
leads to an alternating sequence of -Bfcj's and S*j,'s. As we had seen, the subscripts in this sequence must coincide as 
for the trace of a matrix product. We can also absorb the encounter factor from (|57|) into the latter sequence. The 
indices in that factor represent the original and partner pseudo-orbits to which the first entrance port belongs and 
thus must coincide with the subscripts of the first B in the sequence above. If we sum over all assignments of ports 
to pseudo-orbits by summing over all independent indices, each Z-encounter then gives rise to a factor 

E 



with 



K^n+^'k,)Bk,.nB,,^k,Bk,,,,...Bk„j,B,,,k,^itre{BBf +iXve'{BBf (C7) 
ji,... ,ji,fei,...,fci 



e = diag(ei,e2, . . . ,6^) = diag(eA,ec, • ■ ■ ,ec) 



e 



diag(ei, £2, . . . , e^) = diag(eB, e^,, . . . , ej,) . (C8) 



If we proceed in this way, and also write the exponentials of Eq. (jC6p as J|^.^ e' ''^J'*''^'=''^2+l-S'»3l^) = 

gitr?(i + BS)+itr£'(i+SB)^ ^ ^^^^^ -^^^^ 

Z^r'^= J rf[S,S]e'*""(^+^'^)+'*"~'(5+SB^) J2 ^ J](itre(BB)'(-) +itre'(SB)'M) . (C9) 

^ cnc (7 

Encounter sum: The remaining summation over the number V of encounters and over their sizes /(cr), cr = 
1, 2, . . . , y, is trivial. We get 

ZW(r) = / d[B,B]e'*'-^"(5+stB)+itr.-'(i+SBt)^i. ( (^itr e{BBy +itre'iB By^\ 

V ' \/=2 / 



d[B, 13] exp 




(CIO) 
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2. The sigma model of RMT 



Eq. (jC10|) agrees with the random-matrix result, which we shaU derive in the following. In random- matrix theory 
the replica trick amounts to replacing the spectral determinants in the numerator of the generating function by powers 
A(£')~^''~^) with the number of replicas r an integer larger than unity. One evaluates the thus generalized generating 
function for r > 1, but in the very end takes the limit r ^ 0. This leads to the following representation of the 
generating function 



Z = lim Z(r) 



Zir) 



A E 



A E 



A E 



2'Kp 



-ir-l) 



-(r-l)\ 



E-^ 



(Cll) 



The advantage of the replica trick is that now only inverted spectral determinants show up, which can be represented 
by Gaussian integrals over ordinary complex variables. If we go through the same steps as for the supersymmetric 
sigma model (average over matrices representing the Hamiltonian, Hubbart-Stratonivich transformation, stationary- 
phase approximation) we obtain the following integral representation for Z{r), 



Zir) 



d[Q] exp I -tr 

■2 = „1 V / 



e 
Oe' 



Q 



(C12) 



where we dropped a proportionality factor converging to 1 in the replica limit r — s- 0. In contrast to the supersym- 
metric case, the matrices Q are now complex 2r x 2r matrices. They have to be diagonalizable bt pseudo-unitary 
transformations T which are now defined by 



rA(i)Tt = Ad). 



(C13) 



where A*^^-* = (1,1,..., —1, —1, . . .). The saddle point manifold = —1 has the shape of a hyperboloid. In the replica 
case the relevant solutions of = — 1 are matrices that can be obtained through conjugation from iA^) where , 

Q = TA(^^T-K (C14) 

Multiplication with a block diagonal matrix as in the supersymmetric now leads to the rational parametrization 



T = 1 



B 
B 



(C15) 



where B is an arbitrary r x r matrix and B — B'' . The fact that the off-diagonal blocks of this matrix are Hermitian- 
conjugate to one another follows from the pseudo-unitarity condition. 

In the limit of large e the integral (jC12p can be calculated in a saddle point approximation (not to be mixed with 
the already implemented limit N —^ oo). Using the parametrization (jC15p in Q and expanding the exponent in the 
vicinity of the saddle point _B = we obtain the integral 



Z(i)(r) 



1 d[B,I3]exp 


itre ^ 


1 d[B,i3]exp 





itr." + f^iBBY^ +itri' (^^ + £(55)' j 



1 + BB ,1 + BB 
' e' 



1 - BB 



l-BB 



(C16) 



agreeing with our semiclassical formula. 

Both the semiclassical and the RMT results can be generalized to the orthogonal case. One then obtains a formula 
as m (fCTe]) but with B as a 2r X 2r matrix subject to symmetries analogous to the supersymmetric treatment. In 
fact, even the crossover between the unitary and the orthogonal symmetry class is accessible in this way 44 1. 
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APPENDIX D: FERMIONIC REPLICAS 
1. Semiclassical construction of a sigma model 



A further variant of the sigma model makes use both rephcas and Fermionic variables (the latter appearing only 
in the intermediate stages of the derivation from RMT). To derive this variant from semiclassics, we write the sign 
factor in Eq. ((37|) as 

Here the factor (_ ^^^^^ |-,e dealt with in a similar way as the factor (— ^}jg Bosonic approach. This 
time we split the pseudo-orbits A and B into r — 1 parts. There are (r — 1)"^+"^ ways to do this. If we take the 
limit r ^ this number of possibilities converges to (— 

The remaining sign factor (— i)("B+"o)~("A+nc) depends on the difference between the overall numbers of partner 
orbits and original orbits. To keep track of this factor we envisage the reconnections inside encounters as a step-by- 
step process. In each step the connections of just two stretches are changed, as in Fig. [51 This changes the number 
of orbits by one and thus flips the above sign. If we only have one 2-encounter we need only one such step and the 
sign factor is negative. For the example of a 4-encounter, Fig. [10] shows how we can get from connections in (a) to 
the final connections in (d) in three steps, each time only reshuffling two stretches. In general for each /-encounter 
we need I — 1 steps. Let us denote by vi the number of Z-encounters. Then the overall number of steps, and thus sign 
flips, is given by the number L = '^iVil of encounter stretches or links minus the number V = '^i^i of encounters. 
We are thus free to replace the sign factor by (—1)^^^. 



Altogether we then obtain 



Here the definition of structures is adapted to consider separately situations where the orbits in A or S arc distributed 
differently among the r — 1 subsets of A and B. Furthermore e^, e'j. (and the corresponing diagonal matrices) are now 
defined as 

e = diag(ei,e2, ...,£,.)= diag(eA, ... ,eA,ec) 

e = diag(ei,e2,...,er) =diag(es,---,eB,e£i) • (D3) 

For the diagonal term we have 

= e-^.^^/2-2:.4/2 -Q + (D4) 

j,k—l...r 

If we use (ID2[) as a starting point we can repeat the semiclassical construction of the sigma model from the previous 
appendix. We just have to use the new definitions of e, e' and insert sign factors in the oscillatory prefactor, for each 
of the V encounters, and for each of the L encounter stretch represented by BB, BB. This yields 



Z«(r) = f d[S,S]e-'*""(5-^^)-''"^'(^-^^)^^ l-Y,(itTii-BBy +itTi'i-BBy 



V 



V \ 1=2 



j d[B, B] cxp 



1 ll-BB ,1-BB 
-- tr e = he 

2 \ l + BB l + BB 



(D5) 



2. The sigma model of RMT 

In RMT, the Fermionic replica sigma model is obtained if we represent all spectral determinants the generating 
function ^ by Gaussian Grassmann integrals. To that end one brings the determinants of the denominator "upstairs" 
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FIG. 10: Reconnections in a 4-encounter in three steps 



as [34| det^ = det'^ , with the replica parameter r first taken as an integer larger than unity; only in the end, r 
is assumed real and sent to zero. The spectral determinants can the be represented through Gaussian integrals over 
Fermionic variables. By going through the same steps as in the supersymmetric and the Bosonic replica variant, 
these integrals can be traded for and integral over complex 2r x 2r matrices which are diagonalizable by unitary 
transformations and satisfy = —1, 



Z{r) 



dQe 



itrQ(^_°,) 



(D6) 



/Q2 = _l 

with e, e' as in (|D3p . The matrices Q relevant for the non-oscillatory terms in the correlator are again of the type 

Q = iTK'^^^T-^ . (D7) 

The rational parametrization then amounts to 

-B^ 



T = 1 



B 



which inserted into the integral (|D6p yields 



Z^^\r) = / d[S,S]exp 




1 -BB 
1 + BB 



(D8) 



(D9) 



just as in semiclassics. 

Again both the semiclassical and the RMT calculation can be generalized to the orthogonal case. 

The contribution of the saddle B = B — does not exhaust the high-e asymptotics of the generating function. 
Another contribution is connected with the vicinity of the so called anomalous or Andreev-Altshuler saddle, A^^' = 
diag(— 1, 1, . . . , 1; 1, —1, . . . , —1) which can be obtained from A'^-'^^ by swapping the first and the (r + l)-th entry, i. e., 



A(2) = i?A(i)i?t , 



R 



V 







1 






lr-1 






1 














lr-1 



(DIG) 



where all missing entries are zero. A*-^-* cannot be reached using the rational parametrization (jPSP : its vicinity 
corresponds to B,B tending to infinity. 

To appreciate the contribution of the Andreev-Altshuler saddle, we substitute Q = RQ'R'^ in the integral 
(|D6p and adopt Q' as the new integration variable. The Jacobian of this transformation is unity such that 
we get Zr — jQ,2^idQ' e^*'"'^ ^ an expression still equivalent to (|D6|) . the replacement of X by R'^XR — 



diag(-ei,e2, . 



t2i ■ 



-ej,) notwithstanding. We then employ the rational parametrization for Q' and ob- 



tain the asymptotic series in expanding the integrand around the stationary point B = B = corresponding to 
Q' = A'-'^' or Q = A^^). Again, the Duistermaat-Heckmann theorem guarantees that the series breaks up after the 
leading term. The final result for Z^^'' differs from Z'^^ only by the appearance of the permuted offset matrix R^XR 
instead of X, i. e. 



Z('\eA,eB,ec,eD) = Z(^\ 



(Dll) 



The Weyl symmetry of Z = Z^^^ + Z*^^^ thus arises again. 
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APPENDIX E: STRUCTURES, PERMUTATIONS AND SYMMETRY 

Here we formalize the definition of a structure in terms of permutations, in order to clarify the relation of structures 
to diagrams and investigate the symmetry aspects of that relation. A "diagram" is understood as the set of all 
physically distinct pseudo-orbit quadruplets with a given topology, i.e, with a fixed set of encounters whose ports are 
connected by links in a specified way. Such a diagram can have many equivalent structures associated. For simplicity 
we limit the following to the unitary symmetry class. 

We imagine that the encounters are numbered by cr = 1, . . . , and the stretches and ports in the a-th encounter by 
j — I, . . . , 1{<t). Each port is thus denoted by a label cr.j. According to the convention of section HVBI in the original 
pseudo-orbits A, C each entrance port is connected with the exit port of the same number. After the reconnection 
leading to the pseudo-orbits D the entrance port a.j will be connected to the exit port a.{j — 1), with the exception 
of the first entrance port of each encounter which will be connected to the last exit port. The port reconnections can 
be depicted by a table in which the lower and the upper row list the entrance and exit ports, respectively, and each 
encounter is represented by a block of width ^(cr). 




1.1 ... 1.^(1) 2.1 ... 2.1(2) ... V.l ... V.l{V) \ 
... l.(?(l)-l) 2.1{2) ... 2.(^(2) -1) ... V.l{V) ... V.{l{V)~l)J' 



This table can be viewed as a permutation of L numbers consisting of V independent cycles Pone — 
(1.^(1), . . . , 1.2, 1.1) . . . {V.l{V), . . . V.2, V.l), each cycle associated with an encounter. The port connections by the 
encounter stretches in the original pseudo-orbits A, C are described in this language by the trivial permutation /. 

The connection of exit to entrance ports by links can also be represented by a similar permutation Punk where 
the upper row refers to exit ports and the lower one to entrance ports. Unlike Pcnc the link permutation is totally 
arbitrary. 

The permutation products PunkI = -Piink and Piink-Pcnc define connections between the consecutive entrance ports 
before and after reconnection, respectively. Each cycle of Punk represents the sequence of the entrance ports of a 
pre-reconnection orbit in A, C; similarly, each cycle of -PunkPonc defines such a sequence for a post-reconnection orbit 
in B or D. The numbers n, n' of cycles in these permutations give the numbers of pre- and post-reconnection orbits. 

Our definition of a structure can thus be formalized as a particular choice of (SI) the numbers V and l{a) determining 
the permutation Ponci (S2) a link permutation Punkj and (S3) a division of the n cycles of Punk between the pseudo- 
orbits A and C, and of the n' cycles of PunkPcnc between the pseudo-orbits B and D. 

The order in which the encounters are numbered and the choice of the first stretch in an encounter is irrelevant. 
Their change amounts to port renumbering, i.e., also to a permutation (an arbitrary permutation of the encounter 
numbers, and a cyclic permutation of the labels for each encounter). The possible renumberings corresponding to 
a given V, l{a) form a group C with N{C) = VlYlcr=iK'^) elements. The renumbering C G C transforms the 
encounter and link permutations as P ^ CPC^^ . For a structure devoid of any symmetry, none of the non-trivial 
C G C commute simultaneously with the link and encounter permutation, hence renumberings will produce V\ 11^=1 
structures described by P^'^^ = CPcncC*"-^, Pli^k = C'-PiinkC'"^, all of them corresponding to the same diagram and 
making the same contribution to the generating function. This overcounting is cancelled by the denominator 1/N{C) 
in (|4ip such that we are left with the contribution of the respective diagram in the form of the ratio of the encounter 
and link factors. 

In a symmetric topology C may simultaneously commute with Punk and Pone, and hence the number of equivalent 
structures is smaller than N(C). However renumbering the encounters and stretches will still change the variables 
Sa-jjUuj parametrizing the separations inside the encounters; the indexing of the link durations is also changed. Since 
we integrate over all values of Saj,u^j as well as the link durations, relabelling thus leads to a separate contribution 
even in the symmetric case. To count each quadruplet only once we thus always have to divide by the number of 
possible relabellings N{C). 

Alternatively, we could formulate our approach in terms of a to sum over diagrams rather than structures. However, 
then the diagrams with symmetries would require special treatment. This was done in [i^ for the diagrams with 
L — V = 2, but it becomes too cumbersome for a generalization to arbitrary orders. 

For an example we consider a quadruplet containing, before and after the reconnection, one periodic orbit with two 
2-encounters connected as shown in Fig. 1111 The encounter and link permutations with the port numbering shown in 
Fig. [TT] read 

/l.l 1.2 2.1 2.2\ _ 1-2' 2.1' 2.2'\ 

1^1.2' 1.1' 2.2' 2.1'j ' ""'"1^2.1 2.2 1.2 1.1 j 

where we put primes over the indices of the exit ports, to distinguish them from the entrance ports. We have the 
following possibilities to relabel encounters and stretches (i) interchange all first indices 1 and 2, (ii) interchange 1.1 
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FIG. 11: A symmetric orbit pair with two encounters, already shown in the gallery of Fig. |3]as the uppermost one in the left 
column for L — V = 2; here, one of the possible port numberings is indicated. 



and 1.2, and/or (iii) interchange 2.1 and 2.2. As (i) and simultaneous application of (ii) and (iii) both leave Pcnc 
and Piink invariant there are only two equivalent structures associated with the diagram. However application of 
(i) and simultaneous application of (ii) and (iii) still changes the parameters s^j and Ua-j. Hence when calculating 
the contribution of each structure by integrating over all parameters each quadruplet will be counted four times per 
structure. To avoid avoid overcounting we thus have to divide out N(C) = 8, as we would without any symmetries. 



APPENDIX F: COMPLEX VS REAL CORRELATOR 



Let us prove the relation i?(e) = lim Jj^^^_|^q Re C(e) with the complex correlator defined by in particular the 
presence of the constant 1/2 in Averaging in the definition of the spectral correlator, be it real or complex, is 
carried out over an interval of its both arguments, namely the central energy E and the energy offset e/vrp; to obtain 
a smooth function the product of the two averaging intervals must be large compared with the squared mean level 
spacing. 



5E Se/np:^ {1/pY 



(Fl) 



If we want oscillations of the correlator on the scale of the mean level spacing to be resolved the averaging interval 
of e has to be taken small, <C 1. In view of (IFip this means that the averaging interval [a, b] of the central energy 
with SE = b — a must be large compared with the mean level spacing multiplied by 1/Se. On the other hand the 
averaging interval must be small compared with the range populated by the consecutive levels i?i < . . . < E^. (The 
generalization to unbounded spectra or bounded spectra with iV = oo is straightforward.) We define the average over 

the central energy as (. . .)^ = (6 — a)~^ J^^ dE (...). In the energy range Ei < E < Em the mean level spacing 1/p 
can be taken as constant. Without loss of generality we choose to discuss spectral correlations in the middle of the 
range considered, i. e. at i? = e^+Ei therefore take that energy as the center of the averaging interval [a, b]. 

By explicitly doing the average over E in the definition of the correlator (|4]) we obtain 



N 



^^'^^2 ( 27r2p2 (6 _ a) .£ ll + i?, - 

\ ^ ^ l,k — l TTp * 



-I- -, x—b 

In.""" 



where the remaining average is over e only. Using (x + ir/) ^ ^ V x ^ —iirS (x) for r/ 0^ and lux — In + i Arg a;, 
we can write the real part of C as 



ReC(e) = -i + i?/,/ + i?fl,, 



(F2) 



where Rjj (Rr^r) stem from products of the imaginary (real) parts of the fraction and the square bracket. 

We begin with Ri j assuming Im e+ = ry — > +0. In the summand (i, k) the phase of the expression in the square 
brackets can be (i) 27r (both Ei and E^ lie within the averaging interval of the central energy [a, b]); (ii) tt (only one 
eigenvalue lies within [a, b]); (iii) zero (both eigenvalues are outside [a, b]). Let us assume that the energy offset e/7rp 

is small compared with b — a; then considering the factor S (^-^ + Ei E^j , we may neglect the possibility (ii). Up 
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to an additive constant, Rjj must then coincide with the real correlation function, 



1 



(6 - a)p2 

In the sum of products of real parts 



a<Ei,Ek<b ^ 



E,, 




R{e) + 1. 



Rr,r — 



1 




+ - Ek 



■In 



2-Kp 



E, 



27rp 



Ek 



E^. 



E,, 



(F3) 



the poles of + i?i — i?fe^ are cancelled by the zero of the respective logarithms; the mild logarithmic singularities 
of EiR,R at Ei — a,b ± are unconnected with the level correlation and thus of no physical relevance. (They are 
caused by the rectangular cut-off of the integrand at the borders of the averaging interval in (. . .)^ and disappear 
if this cut-off is smoothed out). Hence the summand has are no relevant singularities or other fluctuations on the 
scale of e. We are thus free to replace the double sum over a discrete eigenvalue spectrum by a double integral over a 
uniformly distributed density, as in 



Rr,r = 



1 



27r2 (6 - a) 



dxdy 



El 



7rp 



In 



Note that the smoothing w.r.t. the offset variable is thus rendered superfluous. 

Consider now the semiclassical limit h ^ oo which means that p ^ cxd; we also assume the spectral span constant (in 
classical units) which means that the number N of levels taken into account goes to infinity. Let the averaging interval 
[a, b] simultaneously shrink compared with the spectral span in such a way that the hierarchy Ej^ — Ei ^ b — a ^ ^/t^P 
is observed. Neglecting ^ and letting the parameter A = ^IZa^ go to infinity we get 



R 



R,R 



A'. 

47r2 A- 



"1 


A 


oo J_ 


-A 


1 dx 






In 


X 





dudv- 



■In 



(l-u) (l + v) 



(l + u) (l-v) 



(F4) 



Collecting (|F2|, jFS]), and (|F4| we come to the desired relation, lim^^o Re C (e+) = i? (e). 
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